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ABSTRACT 

We use giga-particle N-body simulations to study galaxy cluster populations in Hubble volumes of 
ACDM {Vtjn = 0.3, r^A = 0.7) and rCDM {Q,m = 1) world models. Mapping past light-cones of locations 
in the computational space, we create mock sky surveys of dark matter structure to z~ 1.4 over 10, 000 
sq deg and to z~0.5 over two full spheres. Calibrating the Jenkins mass function at z = with samples 
of ~ 1.5 million clusters, we show that the fit describes the sky survey counts to < 20% acccuracy over 
all redshifts for systems more massive than poor galaxy groups (5 x 10^^ /i^^Mo). 

Fitting the observed local temperature function determines the ratio (3 of specific thermal energies 

5/3 

in dark matter and intracluster gas. We derive a scaling with power spectrum normalization (3 oc cTg , 
and find that the ACDM model requires erg = 1-04 to match (3— 1.17 derived from gas dynamic cluster 
simulations. We estimate a 10% overall systematic uncertainty in erg, 4% arising from cosmic variance 
in the local sample and the bulk from uncertainty in the absolute mass scale of clusters. 

Considering distant clusters, the ACDM model matches EMSS and RDCS X-ray-selected survey 
observations under economical assumptions for intracluster gas evolution. Using transformations of mass- 
limited cluster samples that mimic erg variation, we explore Sunyaev-Zel'dovich (SZ) search expectations 
for a 10 sq deg survey complete above 10^** /i~^Mq. Cluster counts are shown to be extremely sensitive 
to (Tg uncertainty while redshift statistics, such as the sample median, are much more stable. Redshift 
information is crucial to extract the full cosmological diagnostic power of SZ cluster surveys. 

For ACDM, the characteristic temperature at fixed sky surface density is a weak function of redshift, 
implying an abundance of hot clusters at z > 1. Assuming constant /3, four kT > 8 keV clusters lie at 
z > 2 and 40 kT > 5 keV clusters lie at z > 3 on the whole sky. Too many such clusters can falsify the 
model; detection of clusters more massive than Coma at z > 1 violates ACDM at 95% confidence if their 
surface density exceeds 0.003 per sq deg, or 120 on the whole sky. 

Subject headings: cosmology:theory — dark matter — gravitation; clusters: general — intergalactic 
medium — cosmology 



1. INTRODUCTION 

Studies of galaxy clusters provide a critical interface 
between cosmological structure formation and the astro- 
physics of galaxy formation. Spatial statistics of the clus- 
ter population provide valuable constraints on cosmolog- 
ical parameters while multi-wavelength studies of cluster 
content offer insights into the cosmic mix of clustered mat- 
ter components and into the interactions between galaxies 
and their local environments. 

In the near future, the size and quality of observed clus- 
ter samples will grow dramatically as surveys in optical, 
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X-ray and sub-mm wavelengths are realized. In the op- 
tical, the ongoing wide-field 2dF (Colless et al. 2001) and 
SDSS (Kepner et al. 1999; Nichol et al. 2001; Annis et 
al. 2001) surveys will map the galaxy and cluster distri- 
butions over large fractions of the sky to moderate depth 
(z ^ 0.3), while deeper surveys are probing of order tens 
of degrees of sky to z~l (Postman et al. 1996; Dalton et 
al. 1997; Zaritsky et al. 1997; Ostrander et al. 1998; Scod- 
deggio et al. 1999; Gal et al. 2000; Gladders & Yee 2000; 
Willick et al. 2001; Gonzalez et al. 2001). In the X-ray, 
ROSAT archival surveys (Scharf et al. 1997; Rosati et al. 
1998; Ebeling et al. 1998; Vikhnin et al. 1998; deCrandi 
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(it al. 1999; Bohringer et al. 2001; Ebcling, Edge & Henry 
2001; Gioia et al. 2001) have generated redshift samples 
of many hundreds of clusters. Similar surveys to come 
from developing Chandra and XMM archives {e.g., Romer 
et al. 2001) will lead to order of magnitude improvements 
in sample size and limiting sensitivity. Finally, the detec- 
tion of clusters via their spectral imprint on the microwave 
background (Sunyaev & Zel'dovich 1972; Birkinshaw 1999) 
offers a new mode of efficiently surveying for very distant 
(z > 1) clusters with hot, intracluster plasma (Barbosa et 
al. 1996; Holder et al. 2000; Kneissl et al. 2001). 

Deciphering the cosmological and astrophysical informa- 
tion in the coming era of large survey data sets requires the 
ability to accurately compute expectations for observables 
within a given cosmology. Given some survey observation 
R at redshift z, a likelihood analysis requires the prob- 
ability p{R, z\C,A) that such data would arise within a 
model described by sets of cosmological C = {Ci} and as- 
trophysical A={Aj} parameters. Considering clusters as 
nearly a one-parameter family ordered by total mass M, 
the likelihood of the observable R can be written 



p{R,z\C,A) 



J dM p{M, z I C) p{R I M, z, A) 
/ dM p{M, z I C) ■ 



(1) 



where p{M, z | C) is the likelihood that a cluster of mass M 
exists at redshift z in cosmology C within the survey of in- 
terest, and p{R \ M, z, A) is the likelihood that observable 
R is associated with such a cluster given the astrophysical 
model A. 

Separating the problem in this way assumes its pieces 
to be independent. The space density n{M, z\C), or mass 
function, describes the probability of finding a cluster at 
redshift z with total mass M within comoving volume el- 
ement dV 

p{M,z\C) oc n{M,z\C)dV. (2) 

The absence of explicit astrophysical dependence in the 
mass function is based on the assiimption that weakly in- 
teracting dark matter dominates the matter energy den- 
sity. If a cluster's total mass M is relatively immune to 
astrophysical processes, then the mass function is well de- 
termined by the gravitational clustering of dark matter. 
On the other hand, the likelihood of a particular observ- 
able feature p{R\M,z,A), is dependent, often critically, 
on the astrophysical model. For optical and X-ray obser- 
vations, it encapsulates the answer to the question "how 
do dark matter potential wells lights up?" 

For Gaussian initial density fluctuation spectra, Press & 
Schechter (1974; PS) used a spherical collapse argument 
and AT-body simulation to show that the space density of 
the rarest clusters is cixponentially sensitive to the am- 
plitude of density perturbations on ~ 10 Mpc scales. 
The analytic form of PS was put on a more rigorous foot- 
ing by Bond et al. (1991), but recent extensions to ellip- 
soidal collapse (Sheth & Tormen 1999; Lee & Shandarin 
1999) revise the original functional form. Calibration by 
N-body simulations has led to a functional shape for the 
mass function that retains the essential character of the 
original PS derivation (Jenkins et al. 2001, hereafter JOl, 
and references therein). For cluster masses defined using 
threshold algorithms tied to the cosmic mean mass den- 
sity Pm{z), JOl show that the mass fraction in collapsed 
objects is well described by a single function that depends 



only on the shape of the filtered power spectrum of initial 
fluctuations (T^(M). 

Complications arise in determining the mass function 
n(M, z) from both simulations and observations. The first 
is semantic. Clusters formed from hierarchical clustering 
do not possess unique, or even distinct, physical bound- 
aries, so it is not obvious what mass to assign to a par- 
ticular cluster. This issue is solvable by convention, and 
we choose here a commonly employed measure Ma, de- 
fined as the mass interior to a sphere within which the 
mean interior density is a fixed multiple A of the critical 
density pdz) at the epoch of identification z. Acknowl- 
edging the non-unique choice of threshold, we develop in 
an appendix a model, based on the mean density profile of 
clusters derived from simulations (Navarro, Frenk & White 
1996; 1997), that transforms the mass function fit param- 
eters to threshold values different from that used here. 

Attempts to empirically constrain the mass function are 
complicated by the inability to directly observe the theo- 
retically defined mass. Instead, a surrogate estimator M 
must be employed that is, in general, a biased and noisy 
representation of Ma. For example, estimates derived 
from the weak gravitational lensing distortions induced on 
background galaxies tend to overestimate Ma by ~ 20%, 
with a dispersion of order unity (Metzler, White & Loken 
2001; White, van Waerbeke & Mackey 2001)). 

The temperature T of the intracluster medium (ICM) 
derived from X-ray spectroscopy is an observationally ac- 
cessible mass estimator. Gas dynamic simulations pre- 
dict that the ICM rarely strays far from virial equilibrium 
(Evrard 1990; Evrard, Metzler & Navarro 1996: Bryan & 
Norman 1998; Yoshikawa, Jing & Suto 2000; Mathiesen 
k Evrard 2001), so that p{M\T) is well described by a 
mean power-law relation with narrow (< 15% in mass) 
intrinsic scatter. Observations are generally supportive of 
this picture (Hjorth, Oukbir & van Kampen 1998; Mohr, 
Mathiesen & Evrard 1999; Horner et al. 1999; Nevalainen 
et al. 2000), but the detailed form of p(M | T) remains un- 
certain. The overall normalization is a particular concern; 
we cannot prove that we know the median mass of, say, a 
6 keV cluster to better than 25% accuracy. 

Even with this degree of uncertainty, the space density 
of clusters as a function of T (the temperature function) 
has been used to place tight constraints on erg, the present, 
linear-evolved amplitude of density fluctuations averaged 
within spheres of radius 8h^^ Mpc. Henry & Arnaud 
(1991) derived as = 0.59 ± 0.02 from temperatures of 25 
clusters in a bright. X-ray flux limited sample, assuming 
rtm = 1- Subsequent analysis of this sample (White, Efs- 
tathiou & Frenk 1993; Eke, Cole & Frenk 1996; Vianna & 
Liddle 1996; Fan, BahcaU & Cen 1997; Kitayama & Suto 
1997; Pen 1998) and revised samples (Markevitch 1998; 
Blanchard et al. 2000) generated largely consistent results 
and extended constraints to arbitrary Qm- For example, 
Pierpaoh, Scott & White (2001), reanalyzing the Marke- 
vitch sample using revised temperatures of White (2000), 
find 



as = 0.495 lo 



0.034 
037 



n 



-0.60 



(3) 



Accurate determination of as is a prerequisite for de- 
riving constraints on the clustered mass density flrn from 
a differential measurement of the local and high redshift 
cluster spatial abundances. Most studies have excluded 
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the possibility that Qm — 1 from current data (Luppino & 
Kaiser 1997; Bahcall, Fan & Cen 1997; Carlberg, Yee & 
Elhngson 1997; Donahue et al. 1998; Eke et al. 1998; Bah- 
call & Fan 1998) but others disagree (Sadat, Blanchard & 
Oukbir 1998; Blanchard & Barlett 1998; Vianna & Liddle 
1999). Uncertainty in trg plays a role in this ambiguity, 
as recently illustrated by Borgani et al. (1999a). In their 
analysis of 16 CNOC clusters at redshifts 0.17 < z < 0.55, 
the estimated value of flm shifts by a factor 3, from 0.35 
to 1.05, as (78 is varied from 0.5 to 0.6. 

Motivated by the need to study systematic effects in 
both local and distant cluster samples, we investigate the 
spatial distribution of clusters in real and redshift space 
samples derived from A^-body simulations of cosmic vol- 
umes comparable in scale to the Hubble Volume (c/Hq)^. 
A pair of 10^ particle realizations of flat cold, dark matter 
(CDM) cosmologies are evolved with particle mass equiv- 
alent to that associated with the extended halos of bright 
galaxies. The simulations are designed to discover the 
rarest and most massive clusters (by maximizing volume) 
while retaining force and mass resolution sufhcient to de- 
termine global quantities (mass, shape, low-order kinemat- 
ics) for objects more massive than poor groups of galaxies 
(~ 5 X IO^^H'^Mq). To facilitate comparison to obser- 
vations, we generate output that traces the dark matter 
structure along the past light-cone of two observing lo- 
cations within the computational volume. These virtual 
sky surveys, along with usual fixed proper time snapshots , 
provide samples of millions of clusters that enable detailed 
statistical studies. We publish the cluster catalogs here as 
electronic tables. 

In this paper, we extend the detailed cluster mass func- 
tion analysis of JOl to the sky survey output, updating 
results using a cluster finding algorithm with improved 
completeness properties for poorly resolved groups. We 
match the observed local X-ray temperature function by 
tuning the proportionality factor /? between the specific 
energies of dark matter and intracluster gas. The required 
value of l3 depends on the assumed as, and we derive a 
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scaling /3~tTg based on virial equilibrium and the Jenk- 
ins' mass function form. From z = subvolumes sized to 
local temperature samples, we show that sample variance 
of temperature-limited samples contributes ^ 4% uncer- 
tainty to determinations of ag- Uncertainty in converting 
temperatures to masses remains the dominant source of 
systematic error in erg, and we investigate the influence of 
a 25 per cent uncertainty in mass scale on expectations for 
Sunyaev-Zel'dovich searches. 

In we describe the simulations, including the pro- 
cess of generating sky survey output, and the model used 
to convert dark matter properties to X-ray observables. 
The cluster mass function is examined in ^|3| Million clus- 
ter samples at 2: = are used to determine the best fit 
parameters of the Jenkins mass function, and we show 
that this function reproduces well the sky survey popula- 
tions extending to z > 1. The interplay between the fit 
parameters, cts and the normalization of cluster masses is 
explored, and this motivates a procedure for transforming 
the discrete cluster sets to mimic variation in erg . 

In §^ we use observations analyzed by Pierpaoli et al. to 
calibrate the specific energy factor (3 for each model. We 
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explore properties of the high redshift cluster population 
in emphasizing uncertainties from erg error, intraclus- 
ter gas evolution and possible X-ray selection biases under 
low signal-to-noise conditions. The dg transformations de- 
veloped in §11 are used to explore cluster yields anticipated 
from upcoming SZ surveys, and the median redshift in 
mass-limited samples is identified as a robust cosmologi- 
cal discriminant. Characteristic properties of the ACDM 
cluster population are summarized in §|^, and we review 
our conclusions in 

2. HUBBLE VOLUME SIMULATIONS 

After an upgrade in 1997 of the Cray T3E at the Rechen- 
zentrum Garching^ to 512 processors and 64Gb of memory, 
we carried out a pair of one bilhon (10^) particle simula- 
tions over the period Oct 1997 to Feb 1999. A memory- 
efficient version of Couchman, Pearce and Thomas' Hydra 
N-body code (Pearce & Couchman 1997) parallehzed us- 
ing shmem message-passing utilities was used to perform 
the computations. MacFarland et al. (1998) provide a de- 
scription and tests of the parallel code. 

We explore two cosmologies with a flat spatial metric, 
a ACDM model dominated by vacuum energy density (a 
non-zero cosmological constant) and a rCDM model dom- 
inated by non-relativistic, cold dark matter. The rCDM 
model completed May 1998 while the ACDM model fin- 
ished Feb 1999. Published work from these simulations in- 
cludes an extensive analysis of counts-in-cells statistics (to 
> lO*'' order) by Colombi et al. (2000) and Szapudi et al. 
(2000), investigation of the clustering behavior of clusters 
(Colberg et al. 2000; Padilla & Baugh 2001), analysis of 
two-point function estimators (Kerscher, Szapudi & Szalay 
2000), a description of the mass function of dark matter 
halos (JOl), a study of confusion on the X-ray sky due to 
galaxy clusters (Voit, Evrard & Bryan 2001), and statistics 
of pencil-beam surveys (Yoshida et al. 2000). Kay, Lid- 
dle & Thomas (2001) use the sky survey catalogs to pre- 
dict Sunyaev-Zel'dovich (SZ) signatures for the planned 
Planck Surveyor mission while Outram et al. (2001) use 
the deep mock ACDM surveys to test analysis procedures 
for the 2dF QSO Redshift Survey. 

2.1. Simulation description 

Table |l| summarizes parameter values for each model, 
including the final epoch matter density fim, vacuum en- 
ergy density r^A, power spectrum normalization erg, start- 
ing redshift Zinu, simulation side length L and particle 
mass m. 

Values of erg were chosen to agree approximately with 
both the amplitude of temperature anisotropics in the cos- 
mic microwave background as measured by COBE and 
with the nearby space density of rich X-ray clusters. The 
degree of uncertainty in these constraints allows the final 
space density of clusters as a function of mass to differ be- 
tween the two simulations. However, as we discuss below, 
it is possible to ensure that the observed space density of 
clusters as a function of X-ray temperature is matched in 
both models by adjusting a free factor [3 used to link X-ray 
temperature to dark matter velocity dispersion. 

To initiate the numerical experiments, particle posi- 
tions and momenta at Zinu are generated by perturbing 
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Table 1 
Model Parameters. 



Model 






0-8 






m}' 


ACDM 


0.3 


0.7 


0.9 


35 


3000 


2.25 


rCDM 


1.0 





0.6 


29 


2000 


2.22 



^Cube side length in h ^ Mpc. 
''Particle mass in 10^^ /i~^Mo. 



a replicated 'glass' of one million particles with a set of 
discrete waves randomly drawn from power spectra com- 
puted for each cosmology. Initial Fourier modes of the ap- 
plied perturbations have amplitudes drawn from a Gaus- 
sian distribution with variance given by the power spec- 
trum P{k) = T{k)Ppiini{k). A Harrison-Zel'dovich pri- 
mordial spectrum Pprim(^0 oc is assumed for both mod- 
els. For the ACDM model, the transfer fimction T{k) is 
computed using CMBFAST (Seljak & Zaldarriaga 1996) 
assuming h = 0.7 and baryon density Oft/i^ = 0.0196 
(Buries & Tytler 1998). The rCDM model uses transfer 
function T{k) = {1 + [aq + {hqf^ + (cg)^] '')-V^, where 
q = k/T, r = Qrnh, a = 6.4 /i"! Mpc, b = Sh'^ Mpc, 
c = 1.7 ft,-^ Mpc and v = 1.13 (Bond & Efstathiou 1984). 

The simulations are designed to resolve the collapse of a 
Coma-sized cluster with 500 particles. Although this res- 
olution is sufficient to capture only the later stages of the 
hierarchical build-up of clusters, convergence tests (Moore 
et al. 1998; Frenk et al. 1999) show that structural prop- 
erties on scales larger than a few times the gravitational 
softening length are essentially converged. From tests pre- 
sented in JOl and in an Appendix to this work, cluster 
identification is robust down to a level of about 20 parti- 
cles. Using Afcoma = l.lxlOiH-iM0 (White et al. 1993), 
leads to particle mass 2.2 x 10^^ /i^^Mq in both mod- 
els, comparable to the total mass within ^ ?,QQh~^ kpc 
of bright galaxies (Fischer et al. 2000; Smith et al. 2001; 
Wilson et al. 2001). The mass associated with one billion 
particles at the mean mass density sets the length L of 
the periodic cube used for the computations, resulting in 
a Hubble Length L = c/Hq = 3/1"^ Gpc for the ACDM 
model and L = 2 /i-i Gpc for tCDM. 

A Newtonian description of gravity is assumed, appro- 
priate for weak-field structures. A non-retarded gravita- 
tional potential is employed because the peculiar acceler- 
ation converges on scales well below the Hubble length. 
The good agreement between the higher-order clustering 
statistics of the simiilations and expectations derived from 
an extended perturbation theory treatment of mildly non- 
linear density fluctuations provides indirect evidence that 
this treatment of gravity is accurate (Colombi et al. 2000; 
Szapudi et al. 2000). 

Gravitational forces on each particle are calculated as 
a sum of a long-range component, determined on a uni- 
form spatial grid of 1024'^ elements using Fast Fourier 
Transforms, and a short-range component found by direct 
summation. The latter force is softened with a spline- 



smoothing roughly equivalent to a Plummer law gravita- 
tional potential (/)(r)cx(r^ + e^)~^/^ with smoothing scale 
e = 0.1/),^-' Mpc. A leapfrog time integration scheme is 
employed with 500 equal time steps for each calciilation. 

Processor time for these computations was minimized 
by employing a parallel algorithm well matched to the 
machine architecture (MacFarland et al. 1998) and by 
simulating large volumes that entail a minimum of mes- 
sage passing overhead. The Cray-T3E offers high in- 
terprocessor communication bandwidth along with a na- 
tive message-passing library {shmem) to control data flow. 
A two-dimensional, block-cyclic domain decomposition 
scheme allocates particles to processors. Each processor 
advances particles lying within a disjoint set of rectangu- 
lar regions of dimension L x {L/32) x {L/16) that subdi- 
vide the computational space. Each calculation required 
approximately 35,000 processor hours, or three days of 
the 512-processor machine. This corresponds to advanc- 
ing roughly 4000 particles per second on an average step. 
The computations were essentially limited by I/O band- 
width rather than cpu speed. Execution was performed in 
roughly twenty stages spanning a calendar time of three to 
four months, with data archived to a mass storage system 
between stages. Approximately 500 Gb of raw data were 
generated by the pair of simulations. 

2.2. Sky survey output 

In addition to the traditional simulation output of snap- 
shots of the particle kinematic state at fixed proper time, 

we introduce her(^ sky survey output that mimics the ac- 
tion of collecting data along the past light-cone of hypo- 
thetical observers located within the simulation volume. 
The method extends to wide-angle surveys an approach pi- 
oneered by Park & Gott (1991) in simulating deep, pencil- 
beam observations. Since there is no preferred location in 
the volume, we chose two survey origins, located at the 
vertex and center of the periodic cube for convenience. 

In a homogeneous world model, a fixed observer at the 
present epoch to receives photons emitted at t < to that 
have traversed comoving distances r{t) = c J^" dta^^{t) 
where a{t) is the scale factor of the metric (a(to) = l) and 
c the speed of light. The set of events lying along the con- 
tinuum of concentric spheres {t, r(t)} for t < to defines the 
past light-cone of that observer. In the discrete environ- 
ment of the numerical simulation, we construct the light- 
cone survey by choosing spherical shells of finite thickness 
such that each particle's state is saved at a pair of consecu- 
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Table 2 
Light-cone Sky Surveys. 



Name 


center 


solid angle 






MS 


(i/2,L/2,L/2) 


in 


0.57 


0.42 


VS 


(0,0,0) 


47r 


0.57 


0.42 


PO 


(0,0,0)- 


n/2 


1.46 


1.25 


NO'' 


{L,L,Lr 


n/2 


1.46 


1.25 


DW 


(0,0,0)- 


10° X 10° 


4.4 


4.6 


XW 


(0,0,0) 


16° X 76° 


6.8 





Orientation centered on cube diagonal. 



''NO/PO have opposite orientations about a common center. 



tive timesteps that bound the exact time of intersection of 
the light-cone with that particle's trajectory. Defining U 
as the proper time at step i of the computation, we choose 
inner and outer radii (1 — r])r(ti^i) and (1 -I- 'r])r{ti_i). 
Here 77 0.02 is a small parameter that safeguards against 
a particle appearing only once in the output record due 
to peculiar motion across the discrete shells during a step. 
The inner radius is set to zero on the final two steps of the 
computation. 

With successive states for particles in the output record, 
a linear interpolation is performed to recover the origi- 
nal second-order time accuracy of the leap-frog integrator. 
Given a particle's position relative to the survey origin Xi 
and, at the subsequent step, a'i+i = Xi + Ax, we solve for 
interpolation parameter a defining position x = Xi + aAxi 
such that |a;| — r{ti + aAt), with At the timestep. For a 
spherical survey, this implies 

a = ""^^-"^ (4) 
2{x{Ax + r{ti)Ar) ^ ' 

with Ar = r(tj) — r(tj-|_i) > 0. After solving for a, the par- 
ticle's position and velocity are interpolated and the result 
stored to create the processed sky survey data sets. Co- 
moving coordinates and physical velocities are stored as 
two-byte integers, sufficient to provide ^ e/10 positional 
accuracy and km/s accuracy in velocity. Data are stored 
in binary form in multiple files, each covering a spatial 
subcube of side L/16 of the entire computational volume. 

The scale L, along with the comoving distance-redshift 
relations shown in Figure 0, determine the redshift extents 
of the surveys listed in Table |^. Two principal survey 
types — spheres and octants — extend to distances L/2 
and L, respectively. From the cube center, the MS full- 
sphere surveys extend to redshifts Zmax = 0.57 (ACDM) 
and 0.42 (rCDM). From the origin and its diagonally op- 
posite image, octant surveys (PO and NO, respectively) 
extend to redshifts 1.46 (ACDM) and 1.25 (rCDM). The 
surveys have opposite orientation; both view the interior 
region of the computational space. The VS sphere centered 
on the origin is created using translational symmetries of 
separate octant surveys conducted from the eight vertices 
of the fundamental cube. The interior portions of the PO 
and NO surveys are thus subsets (opposite caps) of the VS 



survey. The combined volumes of the spheres and octants 
sample the computational volume roughly once for each 
type. In terms of cosmic time, the octants extend over 
the last 74% and 71% of the age of the universe (ACDM 
and rCDM, respectively), equivalent to roughly a 10 Gyr 
look-back time. 

In addition to these surveys, smaller solid angle wedge 
surveys reach to greater depth. A 10 x 10 sq deg deep 
wedge (DW) extends along the cube diagonal to the op- 
posite corner and reaches redshifts 4.4 (ACDM) and 4.6 
(rCDM). For the ACDM model only, a 16 x 76 sq deg 
extended wedge (XW) uses periodic images of the funda- 
mental cube to reach Zmax = 6.8. This wedge is a partial 
extension of the PO survey. 

2.3. Connecting to X-ray observations 

Connecting to observations of clusters requires a model 
that relates luminous properties to the underlying dark 
matter. We focus here on the ICM temperature T under 
the assumption that both T and the dark matter velocity 
dispersion a^, are related to the underlying dark matter 
gravitational potential through the virial theorem (Cava- 
liere & Fusco-Femiano 1976). Empirical support for this 
assumption comes from the observation that T ^ CTgai^ 
(e.g., Wu, Xue & Fang 1999), the scaling expected if both 
galaxies and the ICM are thermally supported within a 
common potential. High resolution simulations of galaxy 
formation within a cluster indicate that cTgai should accu- 
rately reflect the dark matter except for the brightest, 
early-type galaxies which display a mild bias toward lower 
velocity dispersion (Springel et al. 2000). 

Rather than map T to M directly, we prefer to use a 
one-to-one mapping between T and dark matter velocity 
dispersion a^,. This approach has the advantage of natu- 
rally building in scatter between T and M at a fractional 
amplitude ^10% that is consistent with expectations from 
direct, gas dynamic modeling {e.g., Mathiesen & Evrard 
2001; Thomas et al. 2001). Define a one- dimensional ve- 
locity dispersion cr„ by 

^ Np 3 

p i=i j=i 
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where the index i ranges over the Np cluster members iden- 
tified within ta, j sums over principal directions, and v is 
the mean cluster velocity defined by the same Np mem- 
bers. We assume that cr„^ maps directly to X-ray tem- 
perature and introduce as an adjustable parameter the 
ratio of specific energies (3 = a^'^ / {kT / firup) , where k is 
Boltzmann's constant, /i the mean molecular weight of the 
plasma (taken to be 0.59) and nip the proton mass. We fit 
P by requiring that the models match the local tempera- 
ture function. 

Varying the power spectrum normalization erg leads to 
shifts in the space density of clusters as a function of mass 
and velocity dispersion. Values of /? required to fit X-ray 
observations thus depend on erg , an d we derive an approx- 
imate form for this dependence in §4.4 below. 



3. CLUSTER POPULATIONS 

We begin this section by visualizing the evolution of 
clustering in the octant surveys. Details of the cluster 
finding algorithm are then presented, and fits to the z = 
mass functions performed for A = 200. A simple model for 
evolving the fit parameters with redshift in the ACDM case 
is given, and predictions based on the z — fits compared 
to the sky survey output in three broad redshift intervals. 

Additional details are provided in two appendices. Ap- 
pendix A compares output of the SO algorithm employed 
here to that used by JO 4. Appendix B presents a model 
for extending the mass function fits under variation of the 
density threshold A. More generally, it provides a means 
to transform the discrete cluster sample under variation of 
as ■ Appendix C presents the sky survey and z — clus- 
ter catalogs as electronic tables. Truncated versions in the 
print edition list the ten most massive clusters in each sur- 
vey. Electronic versions list all clusters above a mass limit 
of 5 X 40i3/j-i]v[g (22 particles). 

3.4. Evolution of the matter distribution 

In Figure |^, we present maps of the Lagrangian 
smoothed mass density in slices through the octant sur- 
veys that extend to 2; = 4.25. Horizontal and vertical maps 
show comoving and redshift space representations, respec- 
tively. Since the Hubble Length far exceeds the charac- 
teristic clustering length of the mass, the feature most im- 
mediately apparent in the density maps is their overall 
homogeneity. Gravitational enhancement of the cluster- 
ing amplitude over time is evident from the fact that the 
density fluctuations are more pronounced near the survey 
origin (vertex of each triangular slice) compared to the 
edge. The effect is subtle in this image because the dy- 
namic range in density, from black to white, spans three 
orders of magnitude, much larger than the linear growth 
factors of 4.8 (ACDM) and 2.2 (tCDM) for large-scale 
perturbations in the interval shown. 

To verify the accuracy of the clustering evolution in 
the octant surveys, we show in Figure ^ the behavior of 
the rms amplitude of density fluctuations {S^y^^, where 
6 = p/ p„i{z) — 4, within spheres encompassing, on average, 
a mass of 2.2 x 40^^ /i-^M© (4000 particles). Points in the 
flgure show {S'^y^'^ determined by randomly sampling lo- 
cations within twenty radial shells of equal volume in the 
octant surveys. Values are plotted at the volume-weighted 
redshift of each shell. Solid lines are not fits, but show the 
expectations for ((5^)^/^ based on linear evolution of the 



input power spectrum. Deviations between the measured 
values and linear theory, shown in the upper panels of Fig- 
ure ^, are at the 4% level. Although we do not attempt 
here to model these deviations explicitly, Szapudi et al. 
(2000) show that the higher-order clustering properties of 
these simulations at the final epoch are well described by 
an extended perturbation theory treatment of fluctuation 
evolution. 

The orientation of the slice shown in Figure |2| is chosen 
to include the most massive cluster in the ACDM octant 
surveys. It lies at the surprisingly high redshift z = 4.04. 
The inset of Figure ||b shows that this cluster is actively 
forming from mergers fed by surrounding filaments. In 
Figure 0, we show a close-up of the redshift space structure 
in 75, 000 km s~^ wide regions centered aX z = 4.4 and ly- 
ing just interior to the vertical edges of the redshift-space 
views of Figure |^. The grey-scale shows only overdense 
material (5 > and the region includes the most massive 
ACDM cluster. With rest-frame line-of-sight velocity dis- 
persion 4964 km s^^,it produces a ~ 46, 000 km s~^ sized 
'finger-of-God' feature in the lower right image. 

Along with this extreme object, close inspection of Fig- 
ure ^ reveals many more smaller fingers representing less 
massive clusters in the ACDM image. The tCDM panel 
contains many fewer such clusters. It is this difference 
that motivates high redshift cluster counts as a sensitive 
measure of the matter density parameter . To perform 
quantitative analysis of the cluster population, we first de- 
scribe the method used to identify clusters in the particle 
data sets. 



3.2. Cluster finding algorithm 

A number of methods have been developed for identify- 
ing clusters within the particle data sets of cosmological 
simulations. We refer the reader to J04, White (2000) and 
Lacey & Cole (4994) for discussion and intercomparisons 
of common approaches. Two algorithms are employed by 
JO 4. One is a percolation method known as "friends-of- 
friends" (FOF) that identifies a group of particles whose 
members have at least one other group member lying closer 
than some threshold separation. The threshold separation, 
typically expressed as a fraction rj of the mean interparticle 
spacing, is a parameter whose variation leads to families 
of groups, referred to as F0F(7/), with favorable nesting 
properties (Davis et al. 4985). 

The other algorithm of JO 4 is a spherical overdensity 
(SO) method that identifies particles within spheres, cen- 
tered on local density maxima, having radii defined by a 
mean enclosed iso-density condition. We use here an SO 
algorithm that differs slightly from that of JO 4. The iso- 
density condition requires the mean mass density within 
radius ta to be a factor A times the critical density Pc{z) 
at redshift z. JO 4 define spherical regions that are over- 
dense with respect to the mean background, rather than 
critical, mass density. For clarity, we refer to the ap- 
proaches as 'mean SO' and 'critical SO' algorithms. If 
not stated explicitly, reference to SO (A) should be read 
as the critical case evaluated at contrast A. By defini- 
tion, a critical SO (A) population is identical to a mean 
S0(A/f2m) population. 

Our method employs a code independent of that used by 
JO 4. The codes produce matching output for well-resolved 
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Table 3 
Fit parameters for Incr^^(M). 



Model 




a 


b 


ACDM 


0.578 


0.281 


0.0123 


rCDM 


0.527 


0.267 


0.0122 


'^equation 


®- 







groups, but differ at low particle number. Appendix A 
provides a discussion of completeness based on direct com- 
parison of group catalogs from the algorithms. 

3.3. Mass function fits at z^O 

By coadding 22 snapshots of 11 Virgo Consortium 
simulations ranging in scale from 84.5ft.~^Mpc to 
3000 h^^ Mpc, JOl showed that the space density of clus- 
ters defined by either FOF(0.2) or SO(180) algorithms are 
well described by a single functional form when expressed 
in terms of In o'^^(M), where a^{M) is the variance of 
the density field smoothed with a spherical top-hat filter 
enclosing mass M at the mean density. Define the mass 
fraction /(Incr^^) by 



/(Ina-i) 



M dn{< Af , z) 
Pm{z) din (7^1 



(6) 



with pm{z) the background matter density at the epoch of 
interest and n(< M, z) the cumulative number density of 
clusters of mass M or smaller. The general form found by 
JOl for the mass function is 



/(Inrx- 



A exp [- \ \nn-^ +B\^] 



(7) 



where A, B and e are fit parameters. The amplitude A sets 
the overall mass fraction in collapsed objects, plays the 
role of a (linearly evolved) collapse perturbation threshold, 
similar to the parameter Sc in the Press-Schechter model 
or its variants, and e is a stretch parameter that provides 
the correct shape of the mass function at the very dilute 
limit. 

Values of these parameters depend on the particular 
cluster finding scheme implemented, but JOl show them 
to be independent of cosmology and redshift when clus- 
ter masses are based on algorithms tied to the mean mass 
density. For the FOF(0.2) group catalogs, JOl find that 
^ = 0.315, 5 = 0.61 and e = 3.8 provide a fit that describes 
all of the numerical data to < 20% precision over eight 
orders of magnitude in number density. 

We fit here the SO mass function by employing a 
quadratic relation describing the filtered power spectrum 
shape 



In cr-i(M) = -Incris +alnAf + 6(lnAf)^ 



(8) 



where M is mass in units of 10^^/i~^Mo and the rms 
fluctuation amplitude ctis at that mass scale is simply 
related to the fiducial power spectrum normalization by 
IncTis = ln((T8) -|- const. Table ^ lists parameters of the 



fit to equation (g). The maximum error in the fit is 2% 
in Incr"^ for masses above 10^^ h~^MQ. For both models, 
the effective logarithmic slope. 



acff(Af) EE din cr-i(M)/dlnAf 



+ 26 In A/, (9) 



slowly varies between 0.2 and 0.3 for masses in the range 
10^3 to IO^^-^/i-^Mq. The Jenkins mass function (JMF) 
expression for the differential number density as a function 
of mass and redshift 



"-jmf(A/, z) 



M 



-a. 



,ff (Af) exp [ - \lna-\M) + B\ 



(10) 

is the form that we fit to the simulated cluster catalogs. 

The critical SO (200) mass functions at ^ = for both 
models are shown in Figure ^, derived from samples of 
1.39 million (ACDM) and 1.48 million (rCDM) clusters 
above 5 x lO^^/j-^M©. Fits to equation (^) are shown 
as dotted lines, with fit parameters listed in Table ^. The 
upper panels of Figure |g show the fractional deviations 
Sn/n — n/njMF — 1 in bins of width 0.12 in InAf. Er- 
ror bars assume Poisson statistics. For bins with 100 or 
more clusters (Afzoo ^ 2 x IO^^/i-^Mq), the rms devia- 
tions (((5ri/n)2)i/2 are < 3% (Table |). 

The high statistical precision of these fits is a lower 
bound on the absolute accuracy of the mass function cal- 
ibration. Based on the fits performed by JOl to a large 
ensemble of simulations covering a wider dynamic range 
in scale than the HV models alone, we estimate that the 
normalization A may be systematically low by ~ 10% (see 
Appendix A). Considering that this degree of uncertainty 
in A corresponds to an uncertainty in mass of only 2 — 3%, 
this level of accuracy is sufficient for the practical purpose 
of comparing to current and near future observations, for 
which the level of systematic uncertainty in mass is an 
order of magnitude larger. 

Another estimate of systematic error is provided by com- 
paring our results to the recent large simulations of Bode 
et al. (2001). Their 1 h^^ Gpc ACDM simulation assumes 
identical cosmological parameters to our model, but uses a 
particle mass and gravitational softening smaller by factors 
of 30 and 7, respectively. Bode et al. employ a measure of 
mass within a comoving radius 1.5 Mpc; this scale en- 
compasses a critical density contrast of 200 at the present 
epoch for mass M200 = 7.94 X lO^"* /i-^Mq . From their Fig- 
ure 6, the space density of clusters above that mass scale 
is 10"^-^ - 2.5 X lO-^/i^ Mpc"^ In our 27h~^ Gpc^ vol- 
ume, we find 6499 clusters above this mass limit, implying 
a space density 2.4 x I0~'^h^ Mpc"^. 
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Table 4 

so (200) mass function parameters. 



Model 


A 


B 


e 


((5n/n)2)i/2 


ACDM 


0.22 


0.73 


3.86 


0.026 


rCDM 


0.27 


0.65 


3.77 


0.028 



'^equation (p^. 



Figure || shows that the two models do not produce 
identical mass functions at the present epoch; the ACDM 
space density is lower by a factor ^4 than that of tCDM. 
Two factors combine to make this difference. The first is 
that our chosen values of erg straddle the constraint de- 
rived from fitting the local X-ray cluster space density, 
such as that quoted in equation (^). The sense of the 
differences — the ACDM model has lower amplitude and 
rCDM higher, both by about 10% - pushes the models in 
opposite directions. The second factor is that our choice 
of fixed critical threshold A = 200 leads to smaller masses 
for ACDM clusters. Previous work has typically employed 
the lower J7m~dependent thresholds derived by Eke et al. 
(1996) from the spherical collapse solutions of Lahav et al. 
(1991) and Lilje (1992). For f]„ = 0.3. Eke et al. calculate 
critical threshold A = 97.2, leading to masses are larger by 
a factor ^^1.25 compared to A = 200. 

3.4. Sky survey cluster populations 

We define cluster catalogs in the sky survey output us- 
ing the same SO algorithm applied to the snapshots at 
fixed proper time. A minor modification for the ACDM 
model must be made in order to define the threshold A 
with respect to the epoch-dependent critical mass density 

p,(z) = p,(o)[r!™(i + z)3 + r!A]. 

For the choice A = 200, Table |^ lists counts of clus- 
ters identified in the sky survey catalogs above mass lim- 
its 5 X 10", 10^* and IO^/i-^Mq. The lower mass limit 
corresponds to 22 particles and the maximum redshifts of 
the catalogs are given in Table |[ 

Total counts of 1.5 and 0.9 million clusters offer large 
statistical samples. On the other hand, the small numbers 
of objects at the most massive end of the spectrum put the 
finite size of the visible universe into context and provide 
additional motivation for near-future surveys to define an 
absolutely complete sample of the largest clusters in the 
universe. Only a few hundred Coma-like or larger clusters 
are expected on the sky at all redshifts in either cosmology. 

Figure ^ shows redshift-space maps of clusters with 
M200 > lO^'^ft.-^M© and z < 1.25 in 3 x 90 deg strips 
taken from the PO octant surveys of the rCDM (left) 
and ACDM (right) models. The surveys display markedly 
different evolution at z > 0.5; distant clusters are more 
abundant in the low mass density cosmology (Efstathiou, 
Frenk & White 1993; Eke, Cole & Frenk 1996; BahcaU Fan 
& Cen 1997). Within the 3 deg slice — a width equiva- 
lent to two Sloan Digital Sky Survey scans — the ACDM 
model contains 3084 clusters above the lO^*/i~^M0, half 
lying beyond z = 0.70. The rCDM model contains 1122 



clusters, with median 2 = 0.39. 

The sky surface density of clusters within three broad 
redshift intervals are shown as a cumulative function of 
mass in Figure 0. The ranges in redshift are chosen to 
represent three classes of observation: local, z < 0.2; in- 
termediate, 0.2 < z < 0.5, and high, 0.5 < z < 1.2. 
Counts measured within the octant surveys are shown 
as points while solid lines show the number expected 
from integrating the Jenkins mass function N{> M) = 
J^-; dz (dV/dz) X^^, dlnAf'njMF(A/',^). 

For rCDM, the integral is performed using the fit pa- 
rameters determined at z = (Table m. For ACDM, we 
must recognize the fact that, becauseT2m(z) varies along 
the light-cone, the fit parameters will evolve with redshift. 
As rim tends to unity at high redshifts, we expect the 
parameters to converge to the rCDM values. Since differ- 
ences in A, B and e between the two models are small at 
z = 0, we take a simple approach and vary the parameters 
linearly in For example, we assume for A that 



A{nm) = {l-x)A{l) + xA{0.3), 



(11) 



where x={l — r2m)/0.7 and ^(1) and A{0.3) are the z = 
fit parameters for rCDM and ACDM from Table |4[ Simi- 
lar interpolations are assumed for B and e. 

The predictions of this model agree very well with the 
measured counts in the octant surveys. The model is accu- 
rate to < 10% in number for ACDM at all masses and red- 
shifts shown. Similar accuracy is displayed for the rCDM 
model at low and intermediate redshifts, but the model 
systematically underestimates counts in the high redshift 
interval by ^--^25%. 

Dashed lines in Figure |^ show numbers expected by the 
Press-Schechter model in its simplest form (see JOl for de- 
tails). For the rCDM model, the PS curve tends to under- 
estimate the space density at high masses. For the ACDM 
model, the use of mass measured within a critical, rather 
than mean, mass density threshold leads to an offset in 
mass between the measured counts and PS curves at low 
redshifts. The offset declines as ftm approaches unity, re- 
sulting in a relatively good match to the simulated counts 
in the high redshift interval. 

In the ACDM panels, we plot the rCDM JMF curves 
as dotted lines for comparison. At low redshifts, the two 
models exhibit an offset in the direction of rCDM being 
overabundant relative to ACDM, a difference already dis- 
cussed in §3.3 for the z = population. In the intermediate 



redshift interval, this offset is reversed at nearly all masses 
above the 5 x 10^^ h~^MQ hmit. At high redshifts, the 
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Table 5 

counts of clusters with /1m2oo/m0 above mass limits. 







ACDM 






rCDM 




Survey 


5 X 1013 


1014 


1015 


5 X 1013 


1014 


1015 


MS 


564,875 


178,223 


322 


377,043 


102,742 


120 


VS 


565,886 


178,483 


285 


378,548 


103,157 


111 


PO 


255,083 


64,608 


45 


107,900 


22,853 


10 


NO 


259,279 


64,930 


42 


108,807 


23,216 


13 


DW 


5,238 


1,316 


1 


1,833 


411 





Total 


1,504,620 


441,833 


623 


878,356 


226,602 


231 



ACDM counts are typically an order of magnitude higher 
than those of rCDM. 

3.5. Sky survey completeness 

Comparison of the octant counts with the JMF expecta- 
tions provides a measure of the incompleteness of the HV 
sky catalogs that arises from their finite redshift extent. 
Figure ^ plots the cumulative sky surface density of clus- 
ters above mass limits 5 x lOi^, 3 x 10i4, and 10i5 K'^Mq 
as a function of redshift. Points show the sky densities of 
clusters lying at redshift z or higher with masses above the 
stated limits (top to bottom, respectively), determined by 
combining the octant surveys of each model. Lines in the 
figure show the JMF expectations from integrating equa- 
tion ([lo| ) using the linear evolution of the fit parameters, 
equation ([TTI). 

At Coma mass scales (> 10i5 /i^iM©), the catalogs are 
essentially complete, as fewer than one such object is ex- 
pected over TT steradian beyond the survey redshift limit. 
At 3 X 1014 ft,-iM0, the ACDM model octants are missing 
~ 100 clusters expected above z = 1.5, implying 98% 
completeness. The tCDM model at this mass limit is 
essentially complete; the small discrepancy between the 
measured and JMF counts at redshifts z > 0.5 reflects the 
systematic trend exhibited in the high redshift panel of 
Figured. At the mass scale of groups, 5 x lOi^ ft,~iM0, the 
incompleteness becomes more significant. In the ACDM 
model, for example, ~ 15% of the group population should 
lie at z > 1.5. 

The ACDM model possesses a healthy population of 
very high redshift clusters. Across the whole sky, a clus- 
ter as massive as Coma is expected at redshifts as high 
as 1.3. At 2: > 2.5, one 3 x lOi4ft,-iM0 cluster should lie 
somewhere on the sky, joined by ^ 20, 000 others above 
5 X IQi^ h~^MQ, nearly one per square degree. Before get- 
ting carried away by such seemingly firm predictions, we 
must investigate the effect of varying a degree of freedom 
that has so-far been kept fixed in the models: the ampli- 
tude of the fluctuation power spectrum erg. 

4. THE TEMPERATURE FUNCTION, ABSOLUTE MASS 
SCALE AND POWER SPECTRUM NORMALIZATION 
UNCERTAINTY 



In this section, we use the freedom in the mass- 
temperature relation to tune /?, the ratio of specific en- 

^We note a typographical error in their equation (18), which should read ln£ = ^^[(r;i — l)fJ,i + rji ln(l — cxp(— /li))]. 



ergies in dark matter and ICM gas, so that both models 
match the observed local temperature function. We show 
how P should scale with erg so as to maintain consistency 
with observations. The final snapshots are used to cali- 
brate the level of uncertainty in erg arising from sample 
variance in local volumes. We discuss the overall system- 
atic error in erg, then examine in §0 how this uncertainty 
affects predictions for the high redshift cluster population. 

4.1. Fitting local temperature observations 

Pierpaoh, Scott & White (2001, hereafter PSW) pro- 
vide the most recent study of the local temperature func- 
tion and its constraints on erg. The sample of 38 clus- 
ters used in their analysis is adapted from the X-ray 
flux-limited sample of Markevitch (1998) and is designed 
as an essentially volume-limited sample within redshifts 
0.03 < z < 0.09 and galactic latitude |6| > 20° for clusters 
with kT > 6 keV. PSW update temperatures for 23 clus- 
ters in the Markevitch (1998) sample with values given in 
White (2000). 

PSW note that the White temperature values, derived 
from ASCA observations using a multi-phase model of 
cluster cooling flow emission, tend to be hotter (by 14% 
on average), than those Markevitch obtained through a 
single-temperature fit after exclusion of a core emission. 
Based on recent high-resolution studies of cooling flows 
(David et al. 2001; McNamara 2001) that do not appear to 
support the underlying cooling flow emission model used 
by White (2000), there is cause for concern that the in- 
creased temperatures may be artificial. We therefore re- 
vert to the original values of Markevitch (1998), and note 
that the resulting effect on the derived values of f3 is about 
10%. 

We use the data from Tables 3 and 4 of PSW and per- 
form a maximum likelihood fit to determine values of f3 
for each model. Our procedure is similar to that used by 
PSW^, but rather than use an analytic model as a refer- 
ence, we use the binned z = differential velocity distri- 
bution n(a„) converted to a set of temperature functions 
n{T) for /? in the range 0.2 to 2. For a given (3, 300 Monte 
Carlo realizations of the observational sample are gener- 
ated, assuming Gaussian statistics and temperature errors 
distributed evenly in number (half positive, half negative) 
about the best fit value. To consider those clusters for 
which the selection volume is best defined and for which 
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cluster ICM physics is better understood, a lower limit of 
6 keV is applied to each random realization. 

Figure |9| shows cumulative number of clusters as a func- 
tion of temperature. The bold line in each panel gives 
the observations while the thin solid line shows the z — 
snapshot number density obtained using best-fit values 
/3a — 0.92 and /3r = 1.20. The dotted hne shows sky sur- 
vey results for clusters lying within the range observed, 
0.03 < z < 0.09), using combined MS and VS samples and 
values /3a = 0.92 and /3t- = 1.10. 

The fact that a single value of /3a leads to acceptable 
fits for both the snapshot and sky survey samples indicates 
that number evolution in the sky survey sample is small 
for ACDM. The evolution in the tCDM case is sufficient 
to warrant a slightly lower value of Pr for the sky survey 
data. The likelihood analysis described above produces Icr 
error estimates /3a = 0.92±0.06 and ^^ = 1.10±0.06. Cur- 
rent samples constrain the temperature scale of the cluster 
population to an accuracy of about 6%. 

4.2. Cosmic variance uncertainty in erg 

The locally observed cluster sample is one realization of 
a cosmic ensemble that varies due to shot noise and spa- 
tial clustering on survey scales. Except for extending the 
angular coverage of the observations into the plane of the 
galaxy, there is no possibility of observing another cluster 
sample in the same redshift range. Cosmic variance in the 
local sample can be investigated using the large sampling 
volumes of the simulations. 

An impression of the magnitude of the sample varia- 
tion is given in Figure |l^. Differential mass functions 
for 16 independent 5000 sq deg survey volumes, extend- 
ing to z = 0.15 and extracted from the MS and VS sam- 
ples, are shown for the ACDM model. A mass limit of 
5 X 10"'^'* h~^MQ leads to an average sample size of 30 clus- 
ters. Dotted lines in each panel show the range in number 
density of the ACDM Jenkins mass function as as is raised 
and lowered by 5 per cent about its default. The sky sur- 
vey sample variations are largely confined within the ±5% 
range of erg shown. 

The full z = snapshot samples allow a more precise 
estimate of the cosmic variance contribution to erg error. 
We divide the full computational volumes into cubic cells 
of size 375 h-^ Mpc (ACDM) and 400 h-^ Mpc (rCDM). 
Offsetting the grid of cells by half a cell width along the 
principal axes and resampling generates totals of 4096 and 
1000 samples of clusters in cells of volume comparable to 
the 5 X 10^( Mpc)"^ sampled by local temperature ob- 
servations. Within each cell, we determine the most likely 
value of (7g using the JMF space density fit to mass-limited 
samples. Mass limits of 6 x 10^^ (ACDM) and lO^^ h'^MQ 
(rCDM) produce average counts of 30 clusters within each 
cell. 

The distributions of erg resulting from this exercise 
are nearly log-normal; we find erg = 0.911e^''"^° and 
0.591e±° °23 for ACDM and tCDM, respectively. The er- 
ror in erg from the likelihood analysis is well approximated 

by 



Ao-g 
CTg 



((Ar-/V)2)l/2 



N 



erg dn{M) 



n{M) das 



(12) 



above the mass limit, N the mean count, and the Jacobian 
is evaluated at the survey mass limit M^^. The latter is 
a steep function of mass, taking on values 5.8 and 8.3 for 
ACDM and rCDM. 

Because of scatter in the temperature-mass relation, the 
variance of counts-in-cells for temperature- limited samples 
is slightly larger than that of mass-limited samples. Per- 
forming a similar analysis based on counts for temperature 
limited samples results in Aerg/erg = 0.039 for ACDM and 
0.025 for rCDM. Since observations are temperature lim- 
ited, these values apply to analysis of current temperature 
data. 

4.3. M ~ T calibration and overall as uncertainty 

As emphasized by previous studies, uncertainty in the 
calibration of p{M \ T) is the largest source of error in erg. 
The error in as associated with uncertainty in the abso- 
lute mass scale can be derived by solving for the zero in 
the total derivative of the mass function, equation (p^. 
Ignoring the weak mass dependence of Q;cfr(A/), the result 
is 



(^8 



M 



where, at large masses (er(Af) < 



a'{M) 



acfr(M) 



e(ln(T-i(A/) + 5)^-1 



(13) 



(14) 



where {{N - N)^y/^ is the rms deviation of counts-in-cells 



The first term can be connected to a shift in the charac- 
teristic collapsed mass (fixed e^) while the second term, 
which arises from the 1/M factor in equation (Kw, is re- 
quired to maintain constant mass fraction in oBjects at 
fixed lna"^(M). The sensitivity a'{M), plotted in Fig- 
ure |l^, asymptotes to a value 0.4 above ~5 x 10^'*/i~^Mq 
in both cosmologies. Below this mass, a'{M) increases 
considerably, reaching unity at 10^"^ h~^MQ. The rarest, 
most massive clusters place the most sensitive limits on 

erg. 

Attempts at calibrating the mass-temperature relation 
have been made using numerical simulations and observa- 
tions. Simulation results by different groups compiled by 
Henry (2000) and PSW display an overall range of ~ 50% 
in temperature at fixed mass, equivalent to a 75% range 
in mass if one assumes M (x T^/"^. A complicating fac- 
tor is that normalizations are typically quoted using a 
mass-weighted temperature, and this measure can differ 
systematically at the ^ 20% level from the spectral tem- 
peratures derived from plasma emission modeling of the 
simulated ICM (Mathiesen & Evrard 2001). Observational 
attempts at calibrating the relation (Horner, Mushotzky 
& Scharf 1999; Nevalainen, Markevitch & Forman 2000; 
Finoguenov, Reiprich, & Bohringcr 2001) display discrep- 
ancies of similar magnitude to the simulations. Part of 
this variation is due to the fact that these analyses are 
comparing T to estimators M that differ in their degree of 
bias and noise with respect to the theoretical mass Ma. 

With relatively little in the way of firm justification, 
we conservatively estimate the ler uncertainty in the zero 
point of the mass-temperature relation to be AAf/M = 
0.25. Assuming log-normal errors, this assumption allows 
the absolute mass scale to lie within a factor 2.3 range at 
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90% confidence. We note that PSW assume a 15% un- 
certainty in mass, somewhat smaUer than the value used 
here. 

From equation (^3|) in the high mass Umit, the uncer- 
tainty in power spectrum normaUzation is 



0.10 (Icr) 



(15) 



or a 16% uncertainty at 90% confidence. We employ this 
level of error when exploring statistics of high redshift clus- 
ters in 

4.4. Degeneracy in (3 and ag 

The calibration uncertainty discussed above in terms of 
mass can be rephrased in terms of temperature or, equiv- 
alently for this study, the parameter P used to connect 
temperature to dark matter velocity dispersion. An ad- 
vantage of P is that it can be determined independently 
from gas dynamic simulations that model the gravitation- 
ally coupled evolution of the ICM and dark matter. In 
a comparison study of twelve, largely independent simu- 
lation codes applied to the formation of a single cluster, 
Frenk et al. (1999) find good agreement among the com- 
puted specific energy ratios within A = 200, with mean 
and standard deviation /3sim = 1.17 ± 0.05. 

At first glance, this determination agrees well with the 
rCDM value of (3 but is in mild (2.7cr) disagreement with 
the AC DM value derived from the local temperature sam- 
ple in §4.1. However, the uncertainties quoted previously 
for P are derived at the fixed values of as used in the N- 
body simulations. To incorporate the additional sources 
of error in ag discussed above, we use the mass sensitiv- 
ity, equation ([T^), and the virial scaling Tcx/3~^MP with 
p ~ 2/3 exhibited by gas dynamic simulations of clusters 
to derive the scaling 



P oc 



v/a'{M) 



5/3 



(16) 



An increase in P\ sufficient to match the Frenk et al. 
(1999) simulation ensemble value requires erg = 1.04. This 
value is marginally within the range allowed by COBE 
microwave background anisotropy constraints for Hubble 
parameter (Eke, Cole & Frenk 1996). Tighter con- 

straints on (78 could serve to increase the tension between 
the two independent determinations of the specific energy 
ratio. 

5. CLUSTERS AT HIGH REDSHIFT 

We are now in a position to revisit the expected num- 
bers of high redshift clusters, incorporating into the analy- 
sis the systematic uncertainty in power spectrum normal- 
ization. We begin by noting the advantage of predicting 
cluster counts as a function of X-ray temperature rather 
than mass, and compare the model predictions to the sky 
surface density of high redshift clusters from the EMSS 
catalog (Henry et al. 1992; Gioia & Luppino 1994). Red- 
shift information from the RDCS catalog (Rosati et al. 
1998; Borgani et al. 1999b) supports the ACDM model 
under conservative assumptions, but the model predictions 
are sensitive to selection effects related to core luminosity 
evolution. 

We then return to mass selected samples and explore 
the sensitivity of Sunyaev-Zel'dovich searches for distant 



clusters to trg variation. Finally, the redshift evolution of 
characteristic mass and temperature scales at fixed sky 
surface density is used to compare ACDM and rCDM ex- 
pectations against redshift and temperature extremes of 
the observed cluster population. 

5.1. X-ray cluster counts 

Because models are constrained by observations of the 
local temperature function, predictions of counts as a func- 
tion of temperature can be made with smaller uncertainty 
than predictions of counts as a function of mass. The mass 
function requires separate knowledge of /? and erg whereas 
the temperature function requires only a unique combina- 
tion of the pair. This advantage breaks down if P (or an 
equivalent parameter linking T to M) evolves with red- 
shift. Current observations support no evolution (Tran et 
al. 1999; Wu, Xue & Fang 1999), at least for the connection 
between galaxy velocity dispersion and ICM temperature. 
We therefore assume a non-evolving P in order to examine 
the space density of clusters as a function of temperature 
at arbitrary redshift. 

Figure shows the range of cumulative sky counts ex- 
pected as a function of temperature within the same three 
broad redshift intervals used in Figure 0. The range m 
counts shown within each panel corresponds to varying 
In P within its 90% confidence region with erg held fixed for 
each model: 0.83 < Pa < 1.01 and 0.91 < Pr < 1.20. The 
constraint to match local observations produces nearly 
complete overlap in the temperature functions of the two 
cosmologies at z < 0.2. At intermediate redshifts, the 
ACDM counts are boosted by nearly an order of magni- 
tude compared to the low redshift range, while the rCDM 
counts grow by a factor ^ 3. The 90% confidence regions 
for the models become disjoint in this redshift interval. 

In the high redshift region, the models separate further, 
with the characteristic temperature at fixed sky surface 
density a factor ^1.5 times larger in ACDM than rCDM. 
The steep nature of the space density translates this mod- 
erate difference in T into a large factor difference in counts; 
at 8 keV, the counts differ by a factor of about 20. An es- 
timate of the observed sky density in this redshift range, 
based on the EMSS survey data, is shown as the square in 
the upper panel of Figure This point is based on three 
hot {kT > 8 keV) and distant (z > 0.5) clusters covering 
a search area of 278 sq deg (Henry 2000), leading to a sky 
surface density 0.011 per sq deg a,t z > 0.5. The data 
are consistent with the ACDM expectations and rule out 
rCDM at > 99% confidence. 

A modest degree of evolution in Pr with redshift could 
reconcile rCDM with the EMSS data. Additional infor- 
mation, such as the redshift distribution of X-ray flux- 
limited samples, provides independent constraints capa- 
ble of eliminating such a possibility (Oukbir & Blanchard 
1992). The RDCS survey (Rosati et al. 1998) is currently 
the X-ray-selected survey with the most extensive redshift 
data available for distant clusters. The survey, as analyzed 
by Borgani et al. (1999b), is complete within 33 sq deg to 
limiting 0.5 — 2 keV X-ray flux 5 x 10~^^ erg s~^ cm~^ and 
contains 70 clusters with measured redshifts extending to 
near one. 

To explore the compatibility of the octant survey pop- 
ulations with the RDCS sample requires a model for the 
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X-ray luminosity anticipated from the simulated clus- 
ters. As a base model, we assume a mean bolometric 
Lx-T relation = 2.9 x 10^4(T/6 keV)^-^^h-^ erg s'^ 
(Arnaud & Evrard 1999) that is assumed not to evolve 
with redshift (Mushotzky & Scharf 1997; Henry 2000; Fair- 
ley et al. 2000). To account for the fact that the Lx-T 
mapping is not one-to-one, we add a uniformly distributed 
scatter of ±0.4 in \ogn){Lx). Fluxes in an observed 
0.5 — 2 keV X-ray band are derived from a mekal spectral 
synthesis code assuming 0.3 solar metallicity. Applying a 
5 X 10"^"* erg s~^ cm~^ flux cut, excluding z < 0.05 clus- 
ters, and scaling the PO and NO simulated cluster surveys 
to 33 sq deg area leads to predictions shown as the solid 
lines in Figure |l^. Under these economical assumptions of 
ICM evolution, the ACDM model provides an acceptable 
fit to the observations. 

Since the most distant cluster sources are typically de- 
tected at modest signal-to-noise, it is worth investigating 
the influence that additional sources of X-ray emission 
would have on survey selection. The steep nature of the 
mass function offers the opportunity for clusters lying just 
below the survey flux limit to be pushed above it, given 
some mechanism to enhance its X-ray luminosity. 

For the purpose of illustration, we consider adding to the 
base model described above random additional sources of 
X-ray luminosity whose influence increases with redshift. 
These sources may be thought of as arising from cooling 
flows, active galaxies embedded within or near the clus- 
ter, or mergers, any or all of which may be more likely 
at higher redshift. The specific model assumes that half 
of the population has luminosities boosted by an amount 
drawn from a uniform distribution of amplitude a{z)Lx, 
with a{z) — 2z and the base luminosity. Although 
arguably extreme, this model raises the zero-point of the 
Lx-T relation by only 50% at z = 0.5. Expectations for 
RDCS based on this alternative, shown as dashed lines in 
Figure differ significantly from the base model predic- 
tions at high redshift. The rCDM still fails to match the 
observations at redshifts between 0.1 and 0.3, but its high 
redshift behavior is much improved. The ACDM model 
consistently overpredicts the counts beyond z — OA. 

Deep X-ray imaging with Chandra and XMM will help 
settle the issue of whether this toy model is too extreme. 
For now, we note that the good agreement between the 
RDCS and the economical ACDM model predictions may 
signal that the ICM undergoes relatively simple evolution 
dominated by gravitational shock heating after an initial, 
early epoch of preheating (Evrard & Henry 1991; Kaiser 
1991; Bower 1997; Cavahere, Menzi & Tozzi 1999; Balogh, 
Babul & Patton 1999; Llyod-Davies, Ponman & Cannon 
2000; Bower et al. 2001; Bialek, Evrard & Mohr 2001; 
Tozzi & Norman 2001). The preheated cluster simulations 
of Bialek et al. (2001) produce low redshift scaling relations 
for X-ray luminosity, isophotal size and ICM mass versus 
temperature that simultaneously match local observations 
and exhibit little evolution in the Lx-T relation to z^l. 

5.2. Mass-selected samples 

Interferometric SZ surveys have been proposed that 
would survey ~ 10 sq deg of sky per year with sufflcient 
sensitivity to detect all clusters above a total mass limit 

lO^^/i~^M0, nearly independent of redshift (Holder et 
al. 1999; Kneissl et al. 2001). The mass limit assumes that 



the ICM mass fraction does not depend strongly on cluster 
mass or redshift, an assumption supported by simulations. 
Bialek et al. (2001) flnd that the ICM gas fraction within 
A — 200 remains a fair representation of the baryon-to- 
total cosmic ratio: /icm = (0.92 ± 0.04) Qh/ilm above rest 
frame temperature kT = 4 keV. We investigate expecta- 
tions for SZ surveys assuming that they will be sensitive 
to a limiting total mass that is independent of redshift. 

Maps of mass-limited cluster sarnples in SDSS-like sur- 
vey slices were presented in Figure o for the default values 
of (Tg . To illustrate the effect of tJg variation, we plo t clus- 
ters in the same spatial regions again in Figure |lj, after 
applying an effective fractional variation in erg of +10% 
(rCDM) and -10% (ACDM). Ahhough equation (|l|) sug- 
gests a simple shift in mass threshold to mimic a change 
in as, the mass dependence of a'{M) (Figure |ll|) intro- 
duces cumbersome non-linearity into the shift. We adopt 
instead an equivalent procedure that adjusts both masses 
M and number densities n(M) in the HV cluster catalogs 
by amounts 

M' = e^M, 
n(M') d In M' = e"^ n{M) d In M (17) 

with 

\n{l + Aas/as) 

M = 7 ^ U«) 

and (aoff) =0.25. Tests of these transformations using the 
Jenkins mass function verify their accuracy to better than 
10% in number for masses lO^'^'^ — lO^^ "^ h~^MQ and vari- 
ations of power spectrum normalization within the 90% 
confidence region \Aas/crs\ < 0.16. The practical value 
of these simple transformations is in allowing the discrete 
simulation output to represent a family of models covering 
a range of normalizations erg- 

When compared to Figure ^ the intermediate redshift 
cluster populations of the two cosmologies shown in Fig- 
ure ^ appear much more similar. Unlike Figure ^, the 
overall counts above 10^^ h~^MQ in the 3 degree slice are 
now nearly identical — 1696 for rCDM compared to 1843 
for ACDM. However, their redshift distributions remain 
different; the rCDM clusters stay concentrated at lower 
redshifts while the ACDM clusters are more broadly dis- 
tributed (Oukbir & Blanchard 1992). 

Figures || and |lj imply that a redshift statistic, such as 
the sample median, will be superior to counts as a means 
to constrain cosmology. Motivated by the aforementioned 
planned SZ surveys, we perform a specific investigation of 
expectations for a random 10 sq deg survey complete above 
a mass-limit M200 — 10^^/i^^Mq. We sample clusters in 
3000 randomly located, square fields of 10 sq deg area, di- 
vided equally between the PO and NO surveys and chosen 
to avoid survey boundaries. We use the transformations 
in equation dij) to define the cluster population at values 
of CTg different from the default. To drive the models in 
directions that minimize their differences, we increase erg 
in the rCDM model and decrease it in the ACDM case. 

The distributions of counts at all redshifts {z < 1.25), 
counts at high redshift (0.8 <z< 1.25) and the median red- 
shift for clusters above the survey mass limit derived from 
the random samples are presented in Figure |l^. At the de- 
fault values of cg (left column), the distributions of num- 
ber expected either at all redshifts z < 1.25 (bottom row) 
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or at high redshift (middle row) would allow unambiguous 
discrimination between the models using a single 10 sq deg 
field. At high redshift, the ACDM model predicts, on av- 
erage, a factor 15 more clusters than tCDM. Overall, the 
mean counts in 10 sq deg are 117 and 45, respectively. Bi- 
asing as by 10% in the chosen directions (middle column) 
produces essentially identical expectations for the overall 
cluster yield, with both models expecting 72 ± 12 clusters 
per field. At high redshift, the ability to discriminate is 
weakened. For a 16% bias (right column), the sense of the 
overall counts are reversed, with the rCDM model having 
a 60% larger yield, on average, than ACDM. The high red- 
shift count distributions of the models possess considerable 
overlap. 

In contrast to the count behavior, the distributions of 
sample median redshift Zmed are extremely stable to vari- 
ations in as ■ The 95-th percentile value of ^med for rCDM 
moves from 0.498 to 0.528 to 0.538 at 0, 10% and 16% 
bias. As a frequentist measure of discrimination we quote 
the power (Sachs 1982), defined as the probability of re- 
jecting rCDM at the chosen level (95%) of significance 
given ACDM as the true model. Measuring the power 
by integrating the ACDM distributions of z^ncd above the 
95-th percentile rCDM value, results in power of 99.9%, 
98.8% and 94.8%. These power measures, and others cal- 
culated in a similar manner for the counts, arc listed in 
corresponding panels of Figure |l^. High redshift counts 
lose power to discriminate between the models as the ap- 
plied bias on as is increased. 

The large shift in the expected counts as as is varied pro- 
vides an appropriate lever arm to use for placing firmer 
constraints on this parameter with SZ surveys. Holder, 
Haiman and Mohr (2001) estimate that a 10 sq deg survey 
as assumed here could, with complete redshift information 
and assuming complete knowledge of the relation between 
SZ signal and cluster mass, constrain as at the 3 — 5% 
level. 

5.3. Sky surface density of distant clusters 

Chandra X-ray Observatory detections of extended 
X-ray emission from three clusters at z > 1 have re- 
cently been reported. Stanford et al. (2001) report de- 
tection of hot ICM in a pair of RDCS-selected clusters 
separated by only 4 arcmin on the sky and 0.01 in red- 
shift, RX J0848-H4453 at z = 1.27 and RX J0849+4452 
at z = 1.26 (Stanford et al. 1997; Rosati et al. 1999). 
RX J0848-t-4453 appears to have a complex morphol- 
ogy and a cool temperature kT — l.Oto'g while 
RX J0849-)-4452 appears to be a relaxed system with 
higher temperature kT = 5.8^^1'^ keV. In addition to these 
systems, Fabian et al. (2001) present Chandra evidence for 
extended ICM emission at temperature kT = 5.0t keV 
around the radio galaxy 3C294 at z = 1.786. Quoted errors 
in these temperature estimates are 68% confidence values. 

From the temperature-mass relation calibrated by the 
local temperature function sample in §^ and assuming 
a non-evolving /3, we can estimate the masses of these 
clusters. Results for ACDM (rCDM) are 4.4 (3.9) x 
1013/j-iMq and M200 = 3.1 (2.7) x W^'^ h-^MQ for 
RX J0848-h4453 and RX J0849-F4452, respectively, and 
M2oo = l-9 (1.6) X lO^-^/i-^Mg for 3C294. To explore the 
likelihood of finding such clusters, we employ a statistic 
that links physical properties to measurable sky surface 



density. 

The statistics we consider are sky surface density charac- 
teristic mass and temperature, defined as the mass Mjsfi(z) 
and temperature T]s[i{z) at which the differential sky sur- 
face density N'{z) = dN/dz of inversely rank-ordered clus- 
ters at redshift z takes on fixed values. The mass scale 
AIn' {z) is defined by the relation 



N\z) = (1/aurv) 



dlnM n{M,z)dV/dz (19) 



where flsmv is the survey sky area. The characteristic 
temperature is defined in a similar manner. As a practi- 
cal approximation to the redshift differential, we employ 
counts in redshift bins of width 0.1 to derive this statistic 
from the HV sky survey data. 

Figure ^ shows the redshift behavior of the sky surface 
density characteristic (SSDC) mass and temperature for 
the ACDM model. Filled points are values based on the 
combined octant survey populations. Solid lines are pre- 
dictions from the Jenkins mass function, derived by com- 
puting equation ( p^ using equation ([l^) for njMF (M, z) 
and integrating in bins of width 0.1 in redshift. Sky sur- 
face density thresholds dN/dz vary by factors of 10 from 
0.001 to 10 per sq deg per unit redshift, as labeled. Open 
circles show results for the SSDC at 0.01 per sq deg per 
unit redshift extending to z~3 using the 16° x 76° exten- 
sion to the PO survey. Thick dashed lines in each panel 
show the limiting resolved mass of 5 x 10^'^ /i~^Mq (22 
particles) and the corresponding limiting resolved virial 
temperature. The good agreement between the Jenkins 
model and the discrete cluster sample measurements is to 
be expected from the results of Figure 0; the POX exten- 
sion data verify the utility of the model to z~3. 

The vertical bar in each panel of Figure |l^ shows the 
90% uncertainty range in the local calibration of each 
quantity; 6=^° " in kT and 6=^°"^ in M. The HV simu- 
lation and Jenkins model results for the SSDC measures 
can be varied vertically by these amounts in Figure 
The narrow spacings between Mjv/(z) and kTiqi{z) con- 
tours reffect the steepness of the cumulative counts at fixed 
redshift; the terrain of the counts is steep in the mass and 
temperature directions. At a particular redshift, the cali- 
bration uncertainties translate into a large range of allowed 
sky surface densities for a given mass, and a smaller but 
still significant range for a given temperature. 

Although steep in the temperature direction, the con- 
tours in the lower panel of Figure 16 are remarkably flat 
in the redshift direction. Over the entire redshift inter- 
val 0.1 < 2 < 1.5, the JMF expectations for the SSDC 
temperature at 0.01 per sq deg per unit redshift lie in a 
narrow range between 8 and 10 keV. In the ACDM model, 
distant, hot clusters should be as abundant on the sky as 
those nearby. 

Temperatures of the aforementioned observed distant 
clusters are plotted in the lower panel of Figure [l^ as 
open triangles (the RX clusters) and open square (3C294). 
Temperature uncertainties at 90% confidence are shown, 
assuming Gaussian statistics to convert la errors. The 
central values of the hotter pair are consistent with a sky 
surface density of 1 per 10 sq deg per unit redshift, but 
within the temperature measurement uncertainties, these 
objects could be up to a factor 100 more common or a 
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factor^ 1000 more rare. The lower temperature system at 
z — 1.27 is consistent with a surface density of several per 
sq deg per unit redshift. 

Figure ^ shows that the tCDM model is less able to 
accommodate the existence of these z> \ clusters. The 
central temperatures correspond to surface densities of 1 
per 1000 sq deg per unit redshift, a factor 100 times more 
dilute than the ACDM values. Given that only 40 such 
clusters would be expected on the whole sky between red- 
shifts one and two, it would be remarkable that two would 
already be identified by these observations. 

At the most dilute sky surface density plotted in these 
figures, each filled circle represents the hottest or most 
massive cluster within its 0.1-wide redshift bin. Even at 
this highest rank, the variance in the discrete sample SSDC 
values remains remarkably small. An exception is the un- 
usual ACDM a,t z — 1.04. This monster lies nearly a fac- 
tor two above the Jenkins model expectations and its de- 
viation is extreme compared to that displayed by values 
at the same source density and other redshifts. We note 
that its expected temperature of 21 keV exceeds that of 
the hottest known cluster IE 0657-56 at z = 0.296, with 
fcT= 17.4 ± 2.5 keV (Tucker et al. 1998). This cluster is 
the asterisk in Figures |l|and0 

As the hottest known cluster, it is natural to expect 
IE 0657-56 to lie at the extreme end of the surface den- 
sity distribution in the redshift range 0.2 — 0.3. That is 
indeed the outcome of comparing its location to ACDM 
expectations in Figure |l^. For the case of tCDM, its ex- 
istence is more troublesome, but given the combination 
of T calibration uncertainty and scatter demonstrated by 
the first-ranked values of the discrete sample, this system 
is consistent at the ~ 2cr level with the expectations of Fig- 
ure |l^. A similar statement of significance can be made for 
the comparably hot and more distant cluster RX J1347- 
1145, with kT = 14.48li:4^ keV (Ettori, Allen & Fabian 
2001) at z = 0.451. This analysis does not support the 
interpretation of Ettori et al. (2001) that the existence of 
RX J1347-1145 alone can be used to place an upper limit 
on the matter density parameter f2m < 0.5. 

To summarize, interpretation of distant cluster counts 
is complicated by uncertainty in tJg, variation of which can 
lead to large factor changes in yield, as well as uncertainty 
due to possible evolution in j3 and other aspects of astro- 
physical evolution. If a constant (3 assumption is valid for 
ACDM, then --8 keV clusters at z = 1.2 - 1.3 should be as 
numerous on the sky as those lying at z = 0.1 — 0.2. 

6. SUMMARY OF ACDM EXPECTATIONS 

Given the increasing likelihood that the ACDM model 
is an accurate representation of our universe (Pryke et al. 
2001; Netterfield et al. 2001), we provide here a brief sum- 
mary and discussion of the characteristics of its cluster 
population. 

Coma-mass systems. The population of clusters with 
Af2oo in excess of 10^^/i~^Mo is potentially numerous, 
but not overwhelmingly so. With as = 0.9, 400 Coma's 
are expected on the whole sky (Figure ^) , but that num- 
ber ranges between 40 and 2000 as erg is varied within its 
90% confidence limits. The median redshift of this sample 



is expected to be 



0.48, nearly independent of trg- 



Detection of Coma equivalents at z > 1 in excess of .003 
per sq deg (~ 120 across the sky) would rule out ACDM at 



95% confidence. A complete sample of these objects could 
be obtained with an all-sky X-ray imaging survey only 
moderately more sensitive than the ROSAT All-Sky Sur- 
vey (Bohringer et al. 2001). Such a survey would be unique 
in being the first to be absolutely complete, meaning com- 
plete in identifying all members of a class of astrophysical 
objects within the finite volume of our past light-cone. 

Hot X-ray clusters. A characteristic feature of the 
ACDM model is that the hottest clusters populate the 
sky at nearly fixed surface density over a broad redshift 
interval (Figure |l6| ). This implies a testable prediction of 
a nearly flat redshift distribution, within z ~ 0.2 — 1, for 
a temperature-limited sample identified in a fixed angular 
survey area. Within the 10,000 sq deg SDSS area, one 
8 keV cluster is expected to lie at z > 2. 

Clusters at z 3. Looking to higher redshifts, clusters 
with M200 > 10^" /j-^M© and rest frame kT > 4 keV (ap- 
parent kT > 1 keV) should exist at the level of one cluster 
per 100 sq deg per unit redshift under the default trg and 
l3 normalizations (Figure |l6|) . Of order one hundred such 
clusters are to be expected within the SDSS survey area 
in the redshift interval 2.5 — 3.5. Of order ten clusters 
will have rest frame kT > 5 keV and z > 3. The vicin- 
ity of bright quasars may be a natural place to search for 
these systems. Verification of a hot ICM at these redshifts 
would benefit from the large collecting area of the planned 
Constellation-X Observatory. 

Clusters at z < 0.5. The SDSS and 2dF optical surveys 
will provide large numbers of clusters selected in redshift- 
space and extending to redshifts z~0.5. Although these 
samples offer an opportunity to place more sensitive con- 
straints on fjg, a number of systematic effects, such as 
biases in the selection process and the mapping between 
properties measured in redshift space (optical richness or 
velocity dispersion) and underlying cluster mass M, must 
first be carefully calibrated. Such systematic effects can 
be profitably studied by combining semi-analytic models 
of galaxy formation with N-body models of dark mat- 
ter halo evolution {e.g., Springel et al. 2001). An X-ray 
imaging survey to bolomctric flux 10^^^ erg s~^ cm~^, ca- 
pable of identifying all clusters with M2Q0 > 10^'^ h~^MQ 
within z — 0.4 (assuming a non-evolving Lx-T relation), 
would provide the ability to separate truly deep potential 
wells from redshift space superpositions of smaller systems 
(Frenk et al. 1990). 

ICM temperature evolution. In predicting that the red- 
shift distribution of hot clusters at fixed sky surface den- 
sity is flat over observationally accessible redshifts, we have 
implicitly assumed that the X-ray temperature and mass 
follow the virial relation Toe {H{z)M2m)'^^'^ ■ K is impor- 
tant to pursue high resolution imaging and spectroscopy 
of known high redshift clusters with Chandra and XMM 
in order to test whether more complex heating and cool- 
ing processes may be occurring, particularly at high red- 
shift. Such processes would affect attempts to determine 
the geometry of the universe through the X-ray size- 
temperature relation (Mohr & Evrard 1997; Mohr et al. 
2000). 

Precise parameter estimation. The accuracy of con- 
straints on tJg from the cluster temperature function is 
fundamentally limited by the error in normalization of the 
mass-temperature relation of hot clusters, equation (13|). 
One per cent errors on as will require knowing the abso- 
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lute mass scale of clusters to better than 3 per cent. This 
challenging prospect is currently beyond the capabilities 
of direct computational modeling and traditional observa- 
tional approaches, such as mass estimates based on hydro- 
static equilibrium. Weak gravitational lensing, especially 
in the form of "field" lensing (see Mellier & Waerbeke 2001 
for a review), appears the most promising approach; for 
example, Hoekstra et aL(2002) find ag = 0.8i;o;i9 at 95% 
confidence from analysis of relatively bright (limiting mag- 
nitude Rc ^ 24) galaxies in CFHT and CTIO fields cov- 
ering 24 sq deg. Imposing such constraints as priors will 
focus future studies on breaking existing degeneracies be- 
tween dark matter /dark energy densities and astrophysical 
evolution. 

7. CONCLUSIONS 

We present analysis of a pair of giga-particle simulations 
designed to explore the emergence of the galaxy cluster 
population in large cosmic volumes of flat world mod- 
els dominated by matter energy density (rCDM) and a 
cosmological constant (ACDM). Besides shear scale, these 
Hubble Volume simulations are unique in their production 
of sky survey catalogs that map structure of the domi- 
nant dark matter component over large solid angles and 
to depths 2; ~ 1.5 and beyond. Application of a spheri- 
cal overdensity (SO) cluster finding algorithm to the sky 
survey and fixed epoch simulation output results in dis- 
crete samples of millions of clusters above the mass scale 
of galaxy groups (5 x h~^M.Q). These samples form 
the basis of a number of studies; we focus here on precise 
calibration of the mass function and on systematic uncer- 
tainties in cosmological parameter determinations caused 
by imprecise determination of the absolute mass scale of 
clusters. A summary of our principal findings is as follows. 

• We calibrate the SO (200) mass function to the Jenk- 
ins form with resulting statistical precision of better 
than 3% in number for masses between 10^^'^ and 
-|^q15.3 /j-IMq. a preliminary estimate of the overall 
theoretical uncertainty in this calibration is approx- 
imately 20%. 

• We fit the local temperature function under the 

assumption that the disordered kinetic energy in 
dark matter predicts the ICM thermal tempera- 
ture, leading to specific energy ratios /3a = (0.92 ± 
0.06) (crs/O.g)^/^ and (3^ = (1.10 ± 0.06) ((78 /0.6)^/^. 
For the ACDM model, a% = 1.04 is required to match 
the value /3sim = 1-17 preferred by gas dynamic sim- 
ulations of ICM evolution. 

• Based on the Jenkins form for the mass function, we 
derive transformations of the discrete cluster sample 
that mimic variation in erg. Using these transfor- 
mations, we show that the redshift distribution of 
mass-limited samples is a more powerful cosmologi- 
cal diagnostic than cluster counts; the median red- 
shift of clusters more massive than 10"'^^ ^."^M© in a 
single 10 sq deg field of a ACDM cosmology can rule 
out rCDM at a minimum of 95% confidence. 

• The ACDM model, under conservative assumptions 
for intracluster gas evolution, is consistent with 



high redshift cluster samples observed in the X-ray- 
selected EMSS and RDCS surveys. 

• The statistics of sky surface density characteristic 
(SSDC) mass and temperature are introduced to 
more naturally account for observational and the- 
oretical uncertainties in measured physical scales. 
The ACDM model predicts flat redshift behavior in 
the SSDC temperature; a randomly chosen 8 keV 
cluster on the sky is nearly equally likely to lie at 
any redshift in the interval 0.2 to 1.2. 

• With (Tg = 0.9, the ACDM model predicts roughly 
400 Coma-mass (10^'"' /i~^M0) clusters across the sky 
at all redshifts, with the most distant lying just be- 
yond z=l. Pushing CTg to its 95% confidence upper 
limit, the ACDM model could accommodate up to 
120 Coma equivalents on the sky at z>l. 

Larger and deeper cluster samples with accurate deter- 
minations of temperature or mass will lead to improved 
constraints on cosmological and astrophysical parameters. 
The developing 2dF and SDSS surveys will provide large 
numbers of clusters with galaxy velocity dispersion Cgai 
serving as a temperature measure and optical richness A^gai 
serving as a surrogate for mass. Gravitational lensing mass 
estimates will also be possible for co-added ensembles of 
clusters (Sheldon et al. 2001). Extracting cosmological in- 
formation from these data will require likelihoods such as 
p{(jy I fjgai) or p{M I A^gai)- The challenge to the theoretical 
community will be to model these likelihoods at a level of 
precision warranted by the large data sets. Almost cer- 
tainly, the theoretical uncertainty associated with this as- 
pect of the modeling will dominate statistical errors, since 
samples of many thousands, perhaps tens of thousands, 
of groups and clusters will be available in the complete 
2dF and SDSS surveys. By imposing external constraints 
on selected parameters and requiring model consistency 
across independent observables {e.g.,, sub-mm, optical and 
X-ray), constraints on cosmological and astrophysical pa- 
rameters can be derived from a number of alternate paths. 

Valuable complementary information is available at 
X-ray wavelengths. An X-ray imaging survey reaching to 
limiting flux ~ 3 x 10"^'' erg s~^ cm~^ in the 0.5 — 6 keV 
band would be capable of detecting a cluster with 6 keV 
rest frame temperature to z = 1 in either of the cosmolo- 
gies studied here, assuming a non-evolving Lx-T relation. 
The redshift distribution of such a sample would be a pow- 
erful cosmological diagnostic, as long as astrophysical evo- 
lution of the ICM could be sufficiently well constrained. A 
program of deep pointed observations with Chandra and 
XMM (to constrain the astrophysical evolution), coupled 
with a deep X-ray imaging survey covering a significant 
portion of the SDSS area (to identify a large cluster popu- 
lation in redshift space) would be a powerful combination. 
SZ surveys over large solid angle with bolometer arrays 
could play a similar role to an X-ray imaging mission, 
and interferometric arrays will probe to smaller masses 
and higher redshifts than can be achieved by any current 
search techniques. Ultimately, the combination of all these 
approaches, along with deep, optical imaging and spec- 
troscopy, will allow determinations of cosmological param- 
eters to be made not only more precise, by shear statistical 
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weight, but also be made more accurate by improving our 
understanding of the astrophysical processes that govern 
the evolution of the visible components of clusters. 
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See Fig. 2a.gif 



See Fig. 2b.gif 



Fk;. 2. — Maps of the dark matter in slices through the deep octant surveys in (a) rCDM and (b) ACDM world models. Survey origins 
are at the vertices (O), and color represents mass density relative to the mean ranging from 0.05 (black) to 50 (white) on a logarithmic scale. 
Density is mapped onto a two-dimensional grid using an adaptive (Lagragian) smoothing kernel with scale 2 X 10^'^ Ii^^'KIq. For each model, 
two representations of a 45° slice extending to z = 1.25 are shown. Horizontal maps display structure in the comoving metric while the vertical 
maps display the same comoving region, reflected about the diagonal, in redshift space. Positions of clusters at the intersection of filaments 
are evident in redshift space through the radial distortions arising from their internal velocity dispersions (so-called 'fingers of God'). The 
inset of the rCDM image shows the relation between comoving distance r and redshift z over the range mapijcd by the images. The inset of 
the ACDM image shows a close-up of the particle distribution around the largest cluster of the ACDM octant surveys, located aX z = 1.04. 
Particles colored white lie within a sphere of physical radius 1.5 Mpc that encompasses a mean density 200 times the critical value. 
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Fig. 3. — Redshift evolution of the rms amplitude of density fluctuations ((5^)^/^ in top-hat spheres containing, on average, a mass of 
2.2 X 10^^ h~^MQ (1000 particles). Points are octant survey measurements (filled, PO; open, NO) from the ACDM (circles) and tCDM 
(triangles) simulations, obtained by randomly sampling twenty radial shells of equal comoving volume and plotted at the volume- weighted 
redshift of each shell. Solid lines are predictions from linear theory based on the input fluctuation spectra. The upper panel demonstrates 
agreement with linear theory at the ^ 1% level, except for the non-linear departure of rCDM fluctuations at late times. The dashed line in 
the lower panel shows the evolution that ACDM fluctuations would have if they followed the tCDM linear growth evolution. 
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Fig. 4. — Redshift-space structure in 75,000 km wide regions centered at 2: = 1.1 in tCDM (left) and ACDM (right). The grey-scale 
shows only overdense material (5 > 0. The ACDM image includes the most massive cluster in the octant survei-s, visible as the long streak at 
the lower right. The regions shown lie just interior to the vertical edges of the redshift-space maps of Figure eI 
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Fig. 5. — The lower panel shows critical SO(200) mass functions derived frmn the z = HV cluster catalogs (solid lines) along with fits 
to the Jenkins' mass function, equation (lid), using parameters listed in Table H (dotted lines). Upper panels show the percent deviation in 
number density between the HV data and the fits. Error bars are based on Poisson statistics in each mass bin. 
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Fig. 7. — Cumulative sky counts of clusters as a function of mass for low, intermediate and high redshift intervals (bottom to top) for tC 
(left) and ACDM (right). Smooth solid lines in each panel give the expectations from integrating the Jenkins mass function, equation (I 
over the appropriate volumes. Points show counts from the PO (filled circles) and NO (open) octant surveys. Dashed lines are stand 
Press-Schechter estimates. In the ACDM panels, dotted lines display the corresponding JMF expectations for the rCDM cosmology. 
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Fig. 8. — Sky surface density of clusters lying at redshift z or higher in ACDM (filled circles, solid lines) and rCDM (filled triangles, 
dotted lines). Points give numbers derived from the combined octant surveys with masses above M200 = 5 X lO^^ (top), 3 X 10^* (middle) 
and 10^^ ?i~^MQ (bottom). Short vertical lines mark the limiting redshifts of the octant surveys. Lines are expectations at each mass limit 
derived from integrating the Jenkins mass function. 
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Fig. 9. — Thick lines show the local temperature function of Pierpaoli et al. (2001) based largely on data of Markevitch (1998). The HV 
simulation expectations, using best-fit values /3a = 0.92 and /3r = 1.20, are shown from the light-cone (dotted) and ^ = snapshot (solid) 
outputs. The former uses clusters within the combined MS and VS surveys lying in the redshift interval 0.03 < z < 0.09. The latter uses the 
cluster population of the entire computational volume. The vertical dashed line in each panel shows the approximate completeness limit of 
the observations. 
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Fig. 10. — Differential mass functions within 16 independent tt/2 steradiaxi_regions extending to 2 = 0.15, derived from the MS and VS 
surveys of the ACDM model. Dashed lines show JMF expectations, equation (tUj), for erg values varied by ±5% about the input value linearly 
evolved to 2 = 0.1. The volume of the samples is comparable to that of the local observed sample used to constrain eg. 




Fig. 11. — The sensitivity a'(M), equation for the ACDM (soUd) and rCDM (dotted) models. 
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10 

kT [keV] 

Fig. 12. — The range of cumulative sky surface densities of clusters as a function of temperature in three redshift intervals for the ACDM 
(horizontal hatched) and rCDM (vertical) cosmologies. The range is determined from the combined octant survey counts by varying /3 within 
its overall 5 to 95% confidence range. The solid square denotes the EMSS observational result for clusters hotter than 8 keV at 2 > 0.5. 
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Fig. 13. — Points show the redshift distribution of the X— ray flux-Hmited RDCS survey (Rosati et al. 1998; Borgani et al. 1999b), with 
Poissonian errors. SoUd Hues show expectations for a 33 sq deg survey derived from the combined octant surveys by assuming constant /3 and 
a non-evolving Lx—T relation. Dashed lines show plausible confusion effects of core luminosity contamination on the X-ray-selection (see 
text for details). 
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Fig. 15. — The influence of varying erg on the population of clusters more massive than 10^* h~^M0 expected in 10 sq deg survey fields. 
Columns (left to right) show probability distributions at the default erg values, {+/— )10% and (+/— )16% variation (+ for rCDM/— for 
ACDM) for the counts at redshifts z < 1.25 (bottom row), counts in the high redshift interval 0.8 < z < 1.25 (middle) and the median cluster 
redshift (top). Vertical hatched distributions are rCDM, horizontal are ACDM. Numbers in each panel give the power statistic described in 
the text. The distributions are generated by sampling Jij sq deg fields around 3000 randomly chosen pointings in the combined octant surveys 
of each model. Shifts in mass and number, equation (h 7|) are used to effectively vary erg in the cluster catalogs. 
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z 

Fig. 16. — Sky surface density characteristic mass (upper) and temperature (lower) in the ACDM model. Points from the 10,000 sq deg 
combined octant surveys (filled circles) and the 1000 square degree extension (open circles) show values above which the cluster sky surface 
density in the redshift interval z — 0.05 to z + 0.05 exceeds values 0.001, 0.01, 0.1, 1 and 10 per sq deg per unit redshift. Clusters at the 
lowest sky density shown are the most massive or hottest in the particular redshift interval within the combined octants. Note the exceptional 
'monster' cluster at 2 = 1.04. Solid lines are JMF expectations described in the text. Vertical bars in each panel denote the 90% confidence 
uncertainty range in the absolute calibration of mass and temperature scales. Open triangles and squares plot extremes of the known X-ray 
cluster population: RX J0849+4452 at z = 1.26 and RX J0848+4453 at 2 = 1.27 (Stanford et al. 2000, triangles); 3C294 at 2 = 1.786 (Fabian 
et al. 2001, square) and IE 0657-56 at 2 = 0.296 (Tucker et al. 2001, asterisk). Dashed lines show a mass limit of 22 simulation particles 
(upper panel) and the virial temperature at that mass (lower). 
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Fig. 17.— 



Same as Figure hq, but for the rCDM model. 
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APPENDIX 

A. CLUSTER FINDING DETAILS AND COMPLETENESS CHECKS 

The cluster finding algorithm that produced the catalogs in this paper begins by generating a density estimate for each 
particle using the distance to its eighth nearest neighbor (Casertano & Hut 1985), equivalent to Lagrangian filtering on a 
mass scale 2 x 10^^ h~^MQ. Sorting density values in decreasing order provides a list of potential sites for cluster centers. 
The list is pruned by eliminating particles whose densities lie below the threshold Apc{z). Beginning with the first member 
of the sorted list, a sphere of radius ta enclosing mass Ma is defined about that particle so that the enclosed density 
p = 3M/\/4Trr^^ — Apc{z). Particles lying within this sphere are recorded as members of this group and are removed from 
the list of potential cluster centers. The process is repeated sequentially, centering on the next available particle in the 
list ordered by decreasing density, until the list is exhausted. Particles may belong to more than one group, but the center 
of a given group never lies within the spherical boundary of another group. 

In analyzing SO(180) and mean SO(324) populations of the tCDM and ACDM models, respectively, JOl noted a 
problem of incompleteness in the SO cluster finding algorithm at particle counts < 100. Resolution tests in JOl indicated 
that space densities of groups comprised of 20 particles could be underestimated by~30%. We employ here an independent 



SO algorithm with improved completeness properties at small numbers of particles. Figure Al compares the SO(180) and 
mean SO(324) abundance functions (for tCDM and ACDM, respectively) at z = based on the new algorithm to the fits 
published in Appendix B of JOl. 

In the figure, the thick solid line show JOl functional fits while the thin solid and dotted lines show discrete mass functions 
derived with the algorithm employed here. The dot-dashed line is the discrete rCDM mass function derived by JOl using 
the previous SO algorithm. The top panel shows the percent deviation between the discrete sample measurements and 
the fit expectations. 

For the rCDM case, both the old and new algorithms compare well against each other and against the fit above 
~10^^ /i~^Mq. At lower masses, the JOl algorithm displays an increasing underestimate in number density with respect 
to the fit, approaching a 30% underestimate at the mass limit 5 x 10^'^ H^^Mq used in this work. The new SO algorithm 
(dotted) displays a similar qualitative trend, but the underestimate is reduced to < 10% in amplitude. A similar trend is 
seen for the new algorithm in the ACDM case where the number density lies ^^12% lower than the JOl fit expectations. 
This analysis indicates that the amplitude A derived from fitting the space density to the Jenkins form, to equation (0), 
may be biased low by ~10% at masses below ~10^^ h~^MQ. 

A further check of resolution effects is made by directly comparing the HV mass function to one derived from smaller 
volumes with improved mass resolution. We do this for SO(200) clusters in the ACDM model at z — 0, using data from 
the 256"^ particle simulation of a 239.5 Mpc region from Jenkins et al. (1998). The new SO(200) algorithm is used to 
identify clusters in the same manner as done in the HV simulation. The cosmological parameters for the models are the 
same, except for a slight difference in the power spectrum used to generate the initial conditions, whereas the particle 
mass in t he 2 56^ particle simulation is a factor 32 times smaller than that used in the HV computations. 

Figure A2 shows the cumulative number of clusters found in the 256^ particle simulation (dotted line) along with the 



number expected based on the z = HV population (solid). Vertical bars on the HV results show the range in number 
derived from dividing the HV volume into 1728 independent cubes of side 250 Mpc and rescaling the counts in each 
cube to a (239.5 h'^ Mpc)'^ volume. The inset shows the correlation between counts above 5 x 10^^ and 3 x lO^"' H'^Mq 
within the subvolumes. The small- volume simulation result is inconsistent with the HV distribution; the count distribution 
shows agreement at the 98% level at 3 x lO^^/i~^M0 but the 5 x 10^^ count is ~0.5% higher than the maximmn of the 
HV distribution and ~ 20% above the mean. 

These findings, along with the slight discrepancy in predicted versus measured counts in the rCDM octant surveys 
(Figure lead to a conservative estimate of the systematic theoretical uncertainty in the number density of clusters 
above KP^/i^-'^Mq to be 20%. Future studies, in particular those which cross-calibrate results for a particular cosmology 
modeled by different simulation teams, are needed to better assess the overall accuracy of this model of the mass function. 
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Fig. A1. — DifTcrcntial SO mass functions at 2: = are compared to the fits published in JOl. In the lower panel, thick lines are the JOl 
expectations (from their Appendix B) for rCDM at A = 180 (upper curve) and ACDM at mean A = 324 (critical A = 97.2, lower curve). 
The ACDM simulation data are shown by the slightly jagged solid curve. Two simulation results are shown for tCDM — the dot-dashed line 
reflects the SO algorithm used by JOl and the dotted line shows results of the algorithm used in this work. The upper panel displays the 
percent deviation in number density between the HV simulation data and the JOl model fits. 
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Fig. A2. — The cumulative number of clusters within the volume indicated, scaled from the full 2 = HV simulation (solid line) and from 
a single realization of a (239.5 Mpc)-^ volume (Jenkins et al. 1998; dotted line) for the ACDM cosmology. Vertical lines show the entire 
range of counts above masses 5 X 10^"^, 10^* and 3 X 10^^ h~^MQ derived from subsampling 1728 cubic sub-volumes of side 250 Mpc 
within the HV realization and scaling to 239.5 Mpc. The inset plots the correlation of counts above 3 X 10^* (x-axis) and 5 X 10^"^ (y-axis) 
derived from the sub-volumes. The star indicates the Jenkins et al. (1998) values. 
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B. MASS SCALE RENORMALIZATION 



The mass scale of clusters at fixed space density is uncertain, both theoretically and empirically, for reasons discussed in 
the openin g sec tion. The lack of a uniquely defined scale motivates a model that would transform the JMF fit parameters 
derived in §3.3 to values appropriate for a redefined mass scale. As an example, we develop here a mode l to estimate the 
SO mass function fit parameters for threshold values A ^ 200. The method is similar to that used in §^|^ to effectively 
vary ug within the discrete cluster samples. 

We have chosen a convention in which a spherical density threshold A ~ 200 defines cluster masses M. A choice A / 200 
would lead to a new mass M' for each cluster related to the original by some factor 



M' = Me^. 

The space density of a set of disjoint clusters is fixed, implying 

n(Af')dlnM' = n(M)dlnM. 



(Bl) 



(B2) 



This condition, with equation (Bl), when used with the space density, equation (|10[), leads to a relation between JMF 
parameters 



B' 



A' = Aet' 
B-a^a{M)fi 



(B3) 



to first order in fi. 

The factor n can be calculated from an assumed mass profile. We use the form introduced by Navarro et al. (1996; 
hereafter NFW), 



Mix) 



SOOn 



■ Pc r200 



ln(l + cx) — cx/{l + cx) 
ln(l + c) - c/(l + c) 



(B4) 



where x = r/r2oo is a scaled radius and c is a concentration parameter. N-body simulations suggest c~5 at masses near 
10^5 h-^ Mpc for the cosmologies studied here (NFW, Frenk et al. 1999; Jing 2000; Bullock et al. 2001). 

Figure Bl shows the results of applying equation ( p3| ) to critical contrasts A = 97.2 and 500 for the ACDM model, 
assuming c = 5. The logarithmic shifts in mass scale are /i = 0.093 and —0.141, respectively. The agreement between 
the predicted and measured values is quite good. At A = 97.2 (equivalent to the mean contrast of 324 used by JOl), the 
bin-averaged mean fractional error and dispersion (for bins with 10% or less Poisson uncertainty) are only 4.1% and 3.8%. 
At A = 500, the mean is —14.5% and dispersion 4.9%. 
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Fig. B1. — The critical SO(500) and SO(97.2) mass functions (lower and upper, respectively) at 2 = for ACDM. Solid lines are measured 
from the HV simulation while dashed lines are predictions based on rescaling the SO (200) JMF fit (dotted line), assuming an NEW profile 
with concentration parameter c = 5. 
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CLUSTER CATALOGS 

The SO(200) cluster catalogs derived from the sky survey and z — snapshot outputs of the simulations are included 
here as electronic tables. In the print edition, Tables 6 through 17 provide a truncated listing of the ten most massive 
clusters for each survey. Electronic versions list all clusters resolved above a mass limit of 22 particles (5 x 10^^ H'-^Mq), 
counts of which are listed in Table |[ Tables 6 through 11 are ACDM simulation catalogs from the z = snapshot, the 
combined PO and XW sky surveys, and the NO, VS, MS and DW sky surveys, respectively. Tables 12 through 17 provide 
the same for rCDM, with the exception that there is no extended wedge associated with the PO survey of this model (see 
Table I). 

Column entries give the mass M20Q (in 10^^ /i^^M©), redshift (for sky survey data) derived from Hubble flow and the 
radial peculiar velocity, one-dimensional velocity dispersion ay determined from a three-dimensional average (in km s~^), 
position in comoving coordinates (in Gpc), and peculiar velocity (in km s^^). The position is defined as the location 
of the particle having the smallest distance to its eighth neighboring particle (see Appendix A), and the peculiar velocity 
is defined by the mass- weighted mean within r2oo- 
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Table C6 

Clusters in the ACDM z = snapshot survey. 



M200 


a 
km 


X 

/i-i Gpc 


y 

h-^ Gpc 


z 

h-^ Gpc 


km 


Vy 

km 


km 


4.7745 


1795 


2.27792 


0.44894 


0.54195 


-90 


153 


-6 


4.7475 


1866 


2.03410 


0.91683 


0.44397 


269 


46 


727 


4.2840 


1781 


1.74103 


2.84553 


1.40265 


87 


-52 


-154 


3.9645 


1632 


0.96734 


0.65102 


2.32061 


63 


312 


-282 


3.8227 


1891 


1.75213 


0.72151 


1.80426 


380 


-90 


-279 


3.4920 


1665 


1.50312 


1.77336 


2.19206 


-26 


1 


-17 


3.3772 


1531 


1.21292 


2.12486 


0.27272 


-434 


-228 


-821 


3.3592 


1583 


2.50579 


2.18682 


0.31021 


-183 


588 


174 


3.3187 


1519 


2.07928 


2.97238 


1.91753 


-532 


-.531 


-156 


;-!.2(i()2 


15G(i 


0. 38232 


1.70524 


0.62298 


-228 


257 


98 



Table C7 

Clusters in the ACDM PO-l-XW sky survey 



M200 


redshift 


a 


X 


y 


z 


Vx 


V,, 

y 


Vz 






km s^^ 


h-^ Gpc 


h-^ Gpc 


/i-i Gpc 


km s^^ 


km s^^ 


km 


2.1465 


1.04210 


1740 


2.27834 


0.44795 


0.54248 


229 


319 


-99 


1.9035 


0.49167 


1490 


0.48208 


1.20790 


0.20341 


287 


-806 


-725 


1.8337 


0.48731 


1473 


1.25433 


0.30248 


0.12419 


300 


-692 


13 


1.6245 


0.48492 


1440 


1.25607 


0.30110 


0.12310 


-520 


45 


327 


1.6132 


0.74772 


1389 


0.38256 


1.70412 


0.62232 


171 


-58 


192 


1.5052 


0.48217 


1339 


1.00142 


0.47311 


0.64768 


-177 


670 


268 


1.4850 


0.63925 


1309 


1.21259 


0.74437 


0.79747 


-4 


-6 


406 


1.4265 


0.33432 


1338 


0.15125 


0.88516 


0.24984 


-533 


-287 


-210 


1.3927 


0.26386 


1273 


0.62241 


0.37143 


0.16475 


283 


50 


-310 


1.3860 


0.95680 


1403 


1.28228 


0.28511 


1.81830 


-107 


-259 


41 
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Clusters in the ACDM NO sky survey 








M200 


redshift 


(J 


X 


y 


z 


Vx 


Vy 


Vz 






km 


h-^ Gpc 


h-^ Gpc 


h-^ Gpc 


km 


km s ^ 


km s^"'^ 


1.8720 


0.47097 


1357 


0.78598 


0.89484 


0.40376 


-13 


365 


-193 


1.6762 


0.45872 


1332 


0.25973 


1.14905 


0.34105 


-128 


508 


-150 


1.5750 


0.39013 


1393 


0.50888 


0.89409 


0.29059 


-106 


-92 


-379 


1.4985 


0.36176 


1308 


0.54647 


0.73774 


0.40355 


-467 


-203 


-432 


1.4625 


0.09295 


1230 


0.11250 


0.20762 


0.13595 


114 


63 


68 


1.4602 


0.21513 


1254 


0.13535 


0.45678 


0.38800 


-624 


336 


-46 


1.4422 


0.23629 


1209 


0.21785 


0.32614 


0.54597 


41 


299 


-195 


1.4310 


0.84337 


1304 


1.25924 


0.15777 


1.59600 


-158 


-124 


27 


1.3612 


0.25340 


1225 


0.35217 


0.57399 


0.23736 


-38 


458 


-21 


1.35(i7 


0.24869 


1277 


0.()8191 


0.11578 


0.14755 


4(i 


-19 


-1027 
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Table C9 

Clusters in the ACDM VS sky survey 



M200 


redshift 


(7 


X 


y 


z 




Vy 








km s^^ 


h-^ Gpc 


h-^ Gpc 


/i-i Gpc 


km s^^ 


km s^^ 


km s"""^ 


3.7710 


0.36705 


1849 


0.77855 


1.94855 


2.04204 


-179 


109 


-38 


2.3985 


0.12317 


1427 


1.66806 


1.18289 


1.50177 


-100 


-204 


-269 


2.1847 


0.31494 


1513 


1.16682 


1.78942 


0.74591 


-276 


55 


-309 


2.1352 


0.22786 


1510 


1.83355 


0.95563 


1.60533 


228 


-71 


352 


2.0407 


0.25998 


1427 


1.69566 


2.12671 


1.17188 


-529 


84 


-410 


2.0385 


0.27780 


1536 


1.14664 


2.02703 


1.04228 


-73 


94 


220 


1.9957 


0.21552 


1393 


1.20062 


1.88904 


1.87100 


60 


546 


-365 


1.9282 


0.17104 


1473 


1.83406 


1.48791 


1.14040 


367 


-434 


-99 


1.9260 


0.38687 


1515 


1.99907 


1.19847 


2.37855 


294 


-141 


207 


l.!)2.37 


().212;-')2 


142.1 


I.(i2(i8(i 


0.92821 


l.:«r)(i.S 


■m?-, 


-10 


-.IKi 



Table CIO 
Clusters in the ACDM MS sky survey 



M200 


redshift 


(J 


X 


y 


z 




Vy 




lOi5/i~iM0 




km s"""^ 


h-^ Gpc 


h-^ Gpc 


/i-i Gpc 


km 


km s ^ 


km s""'^ 


3.0622 


0.12208 


1765 


1.61387 


1.63047 


1.81743 


-135 


-43 


-404 


2.9452 


0.20045 


1583 


1.13968 


1.94410 


1.59070 


37 


-315 


-414 


2.7225 


0.21132 


1724 


1.77946 


1.14258 


1.88827 


41 


-144 


941 


2.6640 


0.26329 


1622 


1.50191 


1.77218 


2.19369 


434 


317 


-245 


2.3692 


0.39004 


1536 


1.21250 


0.74412 


0.79839 


-210 





684 


2.2320 


0.53759 


1482 


1.20256 


2.12932 


0.26569 


636 


-375 


138 


2.1690 


0.48039 


1534 


1.28226 


0.28481 


1.81857 


55 


-597 


114 


2.0925 


0.54800 


1465 


0.38300 


1.70434 


0.62229 


-215 


29 


239 


2.0272 


0.47799 


1414 


1.16678 


1.69964 


0.28919 


308 


159 


-474 


1.9147 


().4l;-')72 


1;-!!);-! 


0.49259 


1.. --51877 


l.().S747 


8(i 


.148 


-.122 



Table Cll 
Clusters in the ACDM DW sky survey 



M200 


redshift 


(T 


X 


y 


z 


Vx 


Vy 


Vz 






km s^^ 


h-^ Gpc 


h-^ Gpc 


h-^ Gpc 


km 


km s ^ 


km s^"'^ 


1.2667 


0.35165 


1194 


0.52955 


0.51554 


0.62439 


236 


159 


160 


0.7425 


0.63063 


1020 


0.89636 


1.04111 


0.84907 


58 


-57 


21 


0.7335 


0.54100 


1078 


0.83017 


0.88976 


0.74354 


-346 


-269 


-328 


0.6952 


0.86474 


1077 


1.26522 


1.12729 


1.20595 


-417 


-122 


217 


0.6795 


0.23262 


883 


0.42833 


0.36816 


0.34631 


184 


-40 


-201 


0.6277 


0.49817 


968 


0.79520 


0.73876 


0.75309 


24 


340 


-103 


0.6210 


0.61920 


1075 


0.87578 


0.91529 


0.96280 


-130 


-498 


540 


0.5872 


0.16708 


833 


0.29741 


0.28511 


0.25289 


406 


-169 


-255 


0.5557 


0.24148 


766 


0.37255 


0.44021 


0.38042 


331 


-714 


-279 


0.5512 


0.52327 


998 


0.76163 


0.87778 


0.75294 


-389 


-179 


-160 
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Table C12 

Clusters in the tCDM z = snapshot survey. 



M200 


a 
km 


X 

h-^ Gpc 


y 

h-^ Gpc 


z 

hr^ Gpc 


km 


Vy 

km 


km s^^ 


4.8840 


1912 


1.60532 


0.11171 


1.05673 


-273 


241 


-141 


4.8707 


1946 


1.76731 


0.09581 


0.98359 


206 


229 


-728 


4.7597 


1907 


0.46288 


1.83931 


0.73075 


402 


598 


-359 


4.3290 


1787 


1.06007 


1.25750 


0.15836 


195 


-102 


160 


4.3246 


1816 


0.62941 


0.32143 


0.59557 


-120 


462 


-43 


4.1270 


1655 


0.18549 


1.04133 


1.38113 


125 


-584 


-549 


4.0448 


1820 


1.40794 


0.55788 


0.85002 


-596 


282 


-236 


3.8473 


1734 


0.24495 


0.63610 


0.51995 


169 


473 


235 


3.5209 


1688 


1.59185 


0.02998 


1.51849 


-21 


-1044 


-78 


;-').5187 


15(i9 


().547(i9 


().o;-')()i4 


1.24131 


-125 


524 


-()2 



Table C13 
Clusters in the tCDM PO sky survey 



M200 


redshift 


(J 


X 


y 


z 




Vy 


Vz 






km 


h-^ Gpc 


h-^ Gpc 


/i-i Gpc 


km 


km s ^ 


km s""'^ 


2.2333 


0.15351 


1635 


0.02751 


0.34796 


0.21763 


-155 


234 


84 


1.9869 


0.20570 


1568 


0.10704 


0.52032 


0.07637 


-219 


-55 


-49 


1.7671 


0.39740 


1512 


0.62955 


0.31885 


0.59612 


-123 


640 


-105 


1.4208 


0.26994 


1345 


0.18684 


0.32731 


0.55685 


934 


-132 


219 


1.3942 


0.01660 


1259 


0.02224 


0.05106 


0.01156 


-742 


-559 


77 


1.2920 


0.29180 


1311 


0.09007 


0.52090 


0.49083 


-153 


-70 


39 


1.2898 


0.08068 


1246 


0.21347 


0.08099 


0.01491 


122 


-443 


-133 


1.2521 


0.27828 


1341 


0.58237 


0.37616 


0.02971 


-206 


150 


35 


1.1500 


0.25740 


1307 


0.43210 


0.35366 


0.33809 


-507 


-135 


33 


1.0811 


0.19157 


1201 


0.42(i(il 


0.2(i()()() 


0.09128 


-478 


-105 


-KiO 



Table C14 
Clusters in the rCDM NO sky survey 



M200 


redshift 


(T 


X 


y 


z 


Vx 


Vy 


Vz 


lO^^/i-^Mo 




km s^^ 


h-^ Gpc 


h-^ Gpc 


h-^ Gpc 


km 


km s ^ 


km s^^ 


1.8781 


0.16574 


1489 


0.08750 


0.34930 


0.25921 


532 


-14 


-310 


1.6650 


0.23534 


1613 


0.17937 


0.04762 


0.57891 


-102 


-429 


-661 


1.6228 


0.40003 


1540 


0.64304 


0.36559 


0.55891 


-6 


475 


93 


1.5873 


0.23482 


1433 


0.19753 


0.03359 


0.57222 


-426 


-150 


-519 


1.4119 


0.13676 


1456 


0.15879 


0.33140 


0.02784 


70 


401 


403 


1.2388 


0.08092 


1148 


0.13108 


0.17990 


0.05517 


-93 


166 


-475 


1.1566 


0.14827 


1301 


0.23406 


0.24841 


0.20322 


599 


505 


-569 


1.1433 


0.14885 


1309 


0.15585 


0.14787 


0.33972 


36 


-297 


142 


1.1255 


0.13588 


1194 


0.25388 


0.00900 


0.26156 


275 


118 


572 


1.0922 


0.05791 


1091 


0.03287 


0.13919 


0.07741 


23 


353 


197 



Clusters in Hubble Volume Sky Surveys 



Table CI 5 
Clusters in the tCDM VS sky survey 



M200 


redshift 


(7 


X 


y 


z 




y 








km 


h-^ Gpc 


h-^ Gpc 


/i-i Gpc 


km s^^ 


km 


km 


2.5530 


0.08805 


1595 


1.23821 


1.06196 


0.96259 


-18 


-331 


78 


2.5108 


0.15568 


1622 


1.41088 


0.90971 


0.98357 


-261 


-24 


-177 


2.2333 


0.15351 


1635 


1.02751 


1.34795 


1.21762 


-155 


234 


84 


1.9869 


0.20570 


1568 


1.10704 


1.52032 


1.07637 


-218 


-55 


-49 


1.9358 


0.17072 


1389 


0.96303 


1.00973 


0.54585 


-48 


394 


126 


1.8937 


0.16180 


1570 


0.57437 


0.89833 


1.00483 


410 


211 


-184 


1.8781 


0.16574 


1489 


0.91249 


0.65070 


0.74078 


-532 


14 


311 


1.7671 


0.39740 


1513 


1.62955 


1.31885 


1.59612 


-120 


638 


-106 


1.7560 


0.37118 


1547 


1.46201 


0.83718 


1.73270 


127 


852 


-638 


1.7405 


().24.").1() 


1491 


0.92007 


l.l()(i42 


l.r)!)4:-!(i 


-■Mi) 


i;-!? 


81 



Table C16 
Clusters in the tCDM MS sky survey 



M200 


redshift 


(J 


X 


y 


z 




Vy 








km s~^ 


h-^ Gpc 


h-^ Gpc 


/i-i Gpc 


km s^^ 


km s ^ 


km 


2.8616 


0.08778 


1619 


1.12538 


1.09734 


0.80682 


-359 


49 


176 


2.6662 


0.09450 


1569 


0.83169 


1.09487 


0.82472 


-444 


32 


-178 


2.0446 


0.23637 


1552 


0.69721 


1.13044 


1.50301 


161 


-325 


514 


2.0291 


0.23655 


1533 


1.40187 


1.26906 


1.36378 


-248 


-415 


427 


1.9625 


0.36864 


1529 


0.62959 


0.31918 


0.59623 


-93 


539 


-159 


1.8448 


0.37457 


1541 


1.06066 


1.25588 


0.15810 


20 


315 


53 


1.7960 


0.06498 


1418 


1.04302 


1.16623 


0.93822 


566 


268 


62 


1.7804 


0.35602 


1629 


1.35688 


1.63401 


1.44111 


9 


-443 


-117 


1.7649 


0.13530 


1456 


0.70869 


1.17978 


0.85280 


571 


211 


141 


1.7027 


().i;-')i(i7 


i;-!()9 


0.94757 


0.9().S42 


0.65712 


205 


45 


-20 



Table C17 
Clusters in the tCDM DW sky survey 



M200 


redshift 


(T 


X 


y 


z 


Vx 


Vy 


Vz 






km s^^ 


h-^ Gpc 


h-^ Gpc 


h-^ Gpc 


km 


km s ^ 


km s^"'^ 


0.5683 


0.30774 


993 


0.41763 


0.43087 


0.45194 


292 


69 


54 


0.5461 


0.11396 


876 


0.19790 


0.17877 


0.17220 


-94 


-258 


-63 


0.4951 


0.46020 


1001 


0.65398 


0.57027 


0.56948 


-842 


261 


4 


0.4484 


0.32176 


890 


0.47050 


0.45774 


0.42240 


210 


-75 


-31 


0.4285 


0.28315 


837 


0.39105 


0.38532 


0.44474 


37 


-578 


-129 


0.4107 


0.31961 


986 


0.47673 


0.44308 


0.42215 


-146 


375 


11 


0.3929 


0.18746 


818 


0.28105 


0.30862 


0.25462 


351 


600 


-69 


0.3441 


0.17978 


833 


0.28576 


0.26987 


0.26635 


-27 


34 


234 


0.3419 


0.19586 


796 


0.29392 


0.31370 


0.27936 


328 


150 


-391 


0.3397 


0.51005 


894 


0.60047 


0.61408 


0.71184 


70 


119 


190 
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